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We consider chaotic inflation in the theories with the effective potentials which at large <j) behave 
either as <f" or as 6°"^. In such theories inflationary domains containing sufficiently large and 
homogeneous scalar field (j) permanently produce new inflationary domains of a similar type. This 
process may occur at densities considerably smaller than the Planck density. Self-reproduction 
of inflationary domains is responsible for the fundamental stationarity which is present in many 
inflationary models: properties of the parts of the Universe formed in the process of self-reproduction 
do not depend on the time when this process occurs. We call this property of the inflationary 
Universe local stationarity. 

... In addition to it, there may exist either a stationary distribution of probability Pc to flnd a given 

\0 ; field <;/> at a given time at a given point, or a stationary distribution of probability Pp to find a given 

. field at a given time in a given physical volume. If any of these distributions is stationary, we will 

' be speaking of a global stationarity of the infiationary Universe. 

In all realistic infiationary models which are known to us the probability distribution Pc is not 
, stationary. On the other hand, investigation of the probability distribution Pp describing a self- 

^ ' reproducing infiationary Universe shows that the center of this distribution moves towards greater 

' and greater (f> with increasing time. It is argued, however, that the probability of inflation (and of 

the self-reproduction of inflationary domains) becomes strongly suppressed when the energy density 
of the scalar field approaches the Planck density. As a result, the probability distribution Pp rapidly 
approaches a stationary regime, which we have found explicitly for the theories and e"*^. In this 
regime the relative fraction of the physical volume of the Universe in a state with given properties 
(with given values of fields, with a given density of matter, etc.) does not depend on time, both at 
' the stage of infiation and after it. 

, Each of the two types of stationarity mentioned above constitutes a significant deviation of in- 

. fiationary cosmology from the standard Big Bang paradigm. We compare our approach with other 

' approaches to quantum cosmology, and illustrate some of the general conclusions mentioned above 

' with the results of a computer simulation of stochastic processes in the infiationary Universe. 
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I. INTRODUCTION 



The standard Big Bang theory asserts that the Universe was born at some moment t = about 15 billion years ago, 
in a state of infinitely large density and temperature. With the rapid expansion of the Universe the average energy of 
particles, given by the temperature, decreased rapidly, and the Universe became cold. This theory became especially 
popular after the discovery of the microwave background radiation. However, by the end of the 70's it was understood 
that this theory is hardly compatible with the present theory of elementary particles (primordial monopole problem, 
Polonyi fields problem, gravitino problem, domain wall problem) and it has many internal difficulties (flatness problem, 
horizon problem, homogeneity and isotropy problems, etc.). 

Fortunately, all these problems can be solved simultaneously in the context of a relatively simple scenario of the 
Universe evolution — the inflationary Universe scenario 0- The main idea of this scenario is that the Universe 

at the very early stages of its evolution expanded quasi-exponentially (the stage of inflation) in a state with energy 
density dominated by the potential energy density V{4)) of some scalar fleld (j). This rapid expansion made the Universe 
flat, homogeneous and isotropic and decreased exponentially the density of monopoles, gravitinos and domain walls. 
Later, the potential energy density of the scalar field transformed into thermal energy, and still later, the Universe 
was correctly described by the standard hot Universe theory predicting the existence of the microwave background 
radiation. 

The first models of inflation were formulated in the context of the Big Bang theory. Their success in solving internal 
problems of this theory apparently removed the last doubts concerning the Big Bang cosmology. It remained almost 
unnoticed that during the last ten years inflationary theory changed considerably. It has broken an umbilical cord 
connecting it with the old Big Bang theory, and acquired an independent life of its own. For the practical purposes 
of describing the observable part of our Universe one may still speak about the Big Bang, just as one can still use 
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Newtonian gravity theory to describe the Solar system with very high precision. However, if one tries to understand 
the beginning of the Universe, or its end, or its global structure, then some of the notions of the Big Bang theory 
become inadequate. For example, one of the main principles of the Big Bang theory is the homogeneity of the 
Universe. The assertion of homogeneity seemed to be so important that it was called "the cosmological principle" 
0. Indeed, without using this principle one could not prove that the whole Universe appeared at a single moment 
of time, which was associated with the Big Bang. So far, inflation remains the only theory which explains why the 
observable part of the Universe is almost homogeneous. However, almost all versions of inflationary cosmology predict 
that on a much larger scale the Universe should be extremely inhomogeneous, with energy density varying from the 
Planck density to almost zero. Instead of one single Big Bang producing a single-bubble Universe, we are speaking 
now about inflationary bubbles producing new bubbles, producing new bubbles, ad infinitum |3,ll3- Thus, recent 
development of inflationary theory considerably modified our cosmological paradigm • In order to understand better 
this modification, we should remember the main turning points in the evolution of the inflationary theory. 

The first semi- realistic version of inflationary cosmology was suggested by Starobinsky 0]. However, originally it 
was not quite clear what should be the initial state of the Universe in this scenario. Inflation in this model could not 
occur if the Universe was hot from the very beginning. To solve this problem, Zeldovich in 1981 suggested that the 
inflationary Starobinsky Universe was created "from nothing" jllj. This idea, which is very popular now [l^-|l6j|. at 
that time seemed too extravagant, and most cosmologists preferred to study inflation in more traditional context of 
the hot Universe theory. 

One of the most important stages of the development of the inflationary cosmology was related to the old inflationary 
Universe scenario by Guth 3]. This scenario was based on three fundamental propositions: 



1. The Universe initially expands in a state with a very high temperature, which leads to the symmetry restoration 
in the early Universe, (j){T) = 0, where (p is some scalar field driving inflation (the inflaton field). 

2. The effective potential V{(p, T) of the scalar field </> has a deep local minimum at = even at a very low 
temperature T. As a result, the Universe remains in a supercooled vacuum state = (false vacuum) for a long 
time. The energy-momentum tensor of such a state rapidly becomes equal to T^^ — guvViQ), and the Universe 
expands exponentially (inflates) until the false vacuum decays. 

3. The decay of the false vacuum proceeds by forming bubbles containing the field 0o corresponding to the minimum 
of the effective potential V{4)). Reheating of the Universe occurs due to the bubble- wall collisions. 



The main idea of the old inflationary Universe scenario was very simple and attractive, and the role of the old 
inflationary scenario in the development of modern cosmology was extremely important. Unfortunately, as it was 
pointed out by Guth in 'sl, this scenario had a major problem. If the rate of the bubble formation is bigger than the 
speed of the Universe expansion, then the phase transition occurs very rapidly and inflation does not take place. On 
the other hand, if the vacuum decay rate is small, then the Universe after the phase transition becomes unacceptably 
inhomogeneous . 

All attempts to suggest a successful inflationary Universe scenario failed until cosmologists managed to surmount 
a certain psychological barrier and renounce the aforementioned assumptions, while retaining the main idea of ref. Q 
that the Universe has undergone inflation during the early stages of its evolution. The invention of the new inflationary 
Universe scenario ^ marked the departure from the assumptions (2), (3). Later it was shown that the assumption (1) 
also does not hold in all realistic models known so far, for two main reasons. First of all, the time which is necessary 
for the flcld to roll down to the minimum of V{<j), T) is typically too large, so that either inflation occurs before the 
field rolls to the minimum of V{4). T), or it does not occur at all. On the other hand, even if the field (f> occasionally 
was near the minimum of 1^(0, T) from the very beginning, inflation typically starts very late, when thermal energy 
drops from down to V{0,T). In all realistic models of inflation V{0,T) < 10~^°Mp, hence inflation may start 
in a state with = not earlier than at < ~ lO^Mp"^. Duringsuch a time a typical closed Universe would collapse 
before the conditions necessary for inflation could be realized 

The assumption (1) was finally given up with the invention of the chaotic inflation scenario |^. The main idea 
of this scenario was to abandon the assumption that the Universe from the very beginning was hot, and that the 
initial state of the scalar field should correspond to a minimum of its effective potential. Instead of that, one should 
study various initial distributions of the scalar field (p, including those which describe the scalar field outside of its 
equilibrium state, and investigate in which case the inflationary regime may occur. 
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In other words, the main idea of chaotic inflation was to remove all unnecessary restrictions on inflationary models 
inherited from the old theory of the hot Big Bang. In fact, the first step towards this liberation was already made 
when it was suggested to consider quantum creation of inflationary Universe from nothing Chaotic inflation 

naturally incorporates this idea 13], but it is more general: inflation in this scenario can also appear under a more 
traditional assumption of initial cosmological singularity. 

Thus, the main idea of chaotic inflation is very simple and general. One should study all possible initial conditions 
without insisting that the Universe was in a state of thermal equilibrium, and that the field (p was in the minimum 
of its effective potential from the very beginning. However, this idea strongly deviated from the standard lore and 
was psychologically difficult to accept. Therefore, for many years almost all inflationary model builders continued 
the investigation of the new inflationary scenario and calculated high-temperature effective potentials in all theories 
available. It was argued that every inflationary model should satisfy the so-called 'thermal constraints' Tt'I, that 
chaotic inflation requires unnatural initial conditions, etc. 

At present the situation is quite opposite. If anybody ever discusses the possibility that inflation is initiated by 
high-temperature effects, then typically the purpose of such a discussion is to show over and over again that this idea 
does not work (see e.g. [l^), even though some exceptions from this rule might still be possible. On the other hand, 
there are many theories where one can have chaotic inflation (for a review see Q)- 

A particularly simple realization of this scenario can be achieved in the theory of a massive scalar field with the 
effective potential ^f/*^, or in the theory jcj)'^, or in any other theory with an effective potential which grows as (f>" 
at large (jj (whether or not there is a spontaneous symmetry breaking at small cf)) . Therefore a lot of work illustrating 
the basic principles of chaotic infiation was carried out in the context of these simple models. However, it would be 
absolutely incorrect to identify the chaotic inflation scenario with these models, just as it would be incorrect to identify 
new inflation with the Coleman- Weinberg theory. The dividing line between the new inflation and the chaotic inflation 
was not in the choice of a specific class of potentials, but in the abandoning of the idea that the high-temperature 
phase transitions should be a necessary pre-requisite of inflation. 

Indeed, already in the flrst paper where the chaotic inflation was proposed j^, it was emphasized that this scenario 
can be implemented not only in the theories ~ but in any model where the effective potential is sufficiently flat 
at some (j). In the second paper on chaotic inflation this scenario was implemented in a model with an effective 
potential of the same type as those used in the new inflationary scenario. It was explained that the standard scenario 
based on high-temperature phase transitions cannot be realized in this theory, whereas the chaotic inflation can 
occur there, either due to the rolling of the field (j) from > 1, or due to the rolling down from the local maximum 
of the effective potential at (jj = 0. In 20] it was pointed out that chaotic infiation can be implemented in many 
theories including the theories with exponential potentials ~ e""^ with a < ^/T6^^. This is precisely the same class 
of theories which two years later was used in 21] to obtain the power law infiation The class of models where 
chaotic infiation can be realized includes also the models based on the SU{5) grand unified theory 7], the infiation 
(modified Starobinsky model) 23] , supergravity-inspircd models with polynomial and non-polynomial potentials 24j , 
'natural infiation' j25|, 'extended infiation' 'hybrid inflation' |27j] . etc.^ 

Several years ago it was discovered that chaotic inflation in many models including the theories 0" has a very 
interesting property, which will be discussed in the present paper ^]. If the Universe contains at least one 
inflationary domain of a size of horizon (/i-region) with a sufficiently large and homogeneous scalar field (j), then this 
domain will permanently produce new /i-regions of a similar type. During this process the total physical volume of 
the inflationary Universe (which is proportional to the total number of /i-regions) will grow indefinitely. The process 
of self-regroduction of inflationary domains occurs not only in the theories with the effective potentials growing at 
large (j) 9, lOl, but in some theories with the effective potentials used in old, new and extended inflation scenarios as 
well |26. .28.. .291 l30l ]. However, in the models with the potentials growing at large (j) the existence of this effect was 



^ We discussed here this question at some length because recently there appeared many papers comparing observational consequences of 
different versions of inflationary cosmology. Some of these papers use their own definitions of new and chaotic infiation, which differ 
considerably from the original definitions given by one of the present authors at the time when these scenarios were invented. This 
sometimes leads to such misleading statements as the claim that if observational data will show that the inflaton potential is exponential, 
or that it is a pseudo-Goldstone potential used in the 'natural inflation' model, this will disprove chaotic inflation. We should emphasize 
that the generality of chaotic inflation does not diminish in any way the originality and ingenuity of any of its particular realizations 
mentioned above. Observational data should make it possible to chose between models with different effective potentials, such as 
e°"l', the pseudo-Goldstone potential, etc. However, we are unaware of any possibility to obtain inflation in the theories with exponential 
or pseudo-Goldstone potentials outside of the scope of the chaotic inflation scenario. 
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most unexpected, and it leads to especially interesting cosmological consequences |7|. 

The process of the self-reproduction of inflationary Universe represents a major deviation of inflationary cosmology 
from the standard Big Bang theory. In this paper we will make an attempt to relate the existence of this process to 
the long-standing problem of finding a true gravitational vacuum, some kind of a stationary ground state of a system, 
similar to the vacuum state in Minkowski space or to the ground state in quantum statistics. 

A priori it is not quite clear that a stationary ground state in quantum gravity may exist at all. This question 
becomes especially pronounced being applied to quantum cosmology. Indeed, how can one expect that at the quantum 
level there exists any kind of stationary state if all nontrivial classical cosmological solutions are non-stationary? 

A possible (and rather paradoxical) answer to this question is suggested by the observation made in |3lj | that the 
wave function of the whole Universe does not depend on time, since the total Hamiltonian, including the contribution 
from gravitational interactions, identically vanishes. This observation has led several authors to the idea of using 
scale factor of the Universe instead of time, which implied existence of many strange phenomena like time reversal 
and resurrection of the dead at the moment of maximal expansion of a closed Universe. The resolution of the paradox 
was suggested by DeWitt in [sJl (see also discussion of this question in [s^l, H)- We do not ask why the Universe 
evolves in time measured by some nonexistent observer outside our Universe; we just ask why we see it evolving in 
time. At the moment when we make our observations, the Universe is divided into two pieces: an observer with its 
measuring devices and the rest of the Universe. The wave function of the rest of the Universe does depend on the 
time shown by the clock of the observer. Thus, at the moment we start observing the Universe it ceases being static 
and appears to us time-dependent. 

From this discussion it follows that the condition of time-independence is not strong enough to pick up a unique 
wave function of the Universe corresponding to its ground state: each wave function satisfying the Wheeler-DeWitt 
equation (being interpreted as a Schrodinger equation for the wave function of the Universe) is time-independent. 
One way to deal with this problem is to assume that the standard Euclidean methods, which help to find the wave 
function of the vacuum state in ordinary quantum theory of matter fields, will work for the wave function of the 
Universe as well. This assumption was made by Hartle and Hawking 33]. Another way is to look for a possibility 
that with an account taken of quantum effects our Universe in some cases may approach a stationary state, which 
could be called the ground state. 

In some cases these two approaches lead to the same results. For example, the square of the Hartle-Hawking wave 
function correctly describes Pc{(t>,t), the stationary probability distribution for finding the scalar field (jj at a given 
time in a given point of de Sitter space with the Hubble constant H 3> m, where m is the mass of the scalar field 
|34|-|36||. Here the subscript c in Pc{(l),t) means that the distribution is considered in comoving coordinates, which 
do not take into account the exponential growth of physical volume of the Universe. However, later it was realized 
that no stationary solutions for Pc{(t>,t) can exist in realistic versions of inflationary cosmology (E,.10]. The reason is 
that the condition H ^ m is strongly violated near the minimum of the effective potential V{(j)) corresponding to the 
present state of the Universe. Indeed, this condition could be satisfied at present only if the Compton wavelength of 
the inflaton field were larger than the size of the observable part of the Universe ^ H^^ and its mass were smaller 
than 10^^^ eV. In all realistic models of infiation this condition is violated by more than thirty orders of magnitude. 

Fortunately, another kind of stationarity may exist in many models of the inflationary Universe due to the process of 
self- reproduction of the Universe 0,0) 01- The properties of inflationary domains formed during the self-reproduction 
of the Universe do not depend on the moment of time at which each such domain is formed; they depend only on 
the value of the scalar flelds inside each domain, on the average density of matter in this domain and on the physical 
length scale. For example, all domains of our Universe with energy density po lO'^^g • cm^^ filled with the same 
scalar fields look alike on the same length scale, independently of the time when they were formed. This kind of 
stationarity, as opposed to the stationarity of the distribution Pc{(j),t) corresponding to the Hartle-Hawking wave 
function, cannot be described in the minisuperspace approach. Hopefully, the existence of this stationarity (which 
implies that our Universe has a fractal structure) will eventually help us to find the most adequate wave function 
describing the self-reproducing inflationary Universe. 

One way to describe this kind of stationarity is to use the methods developed in the theory of fractals. This question 
was studied in j|2j in application to the theories where inflation occurs near a local maximum of the effective potential 
V{(j)). In this case the expansion of the whole Universe, though eternal, is almost uniform — the scale factor of the 
Universe a{t) grows as e^', where H almost does not depend on the value of the inflaton field 0. This makes it possible 
to factor out the trivial overall expansion factor, and come to a time-independent fractal structure without taking into 
account the difference between Pc{(j),t) and Pp{(j),t). Unfortunately, it is rather difficult to apply these methods to 
the most interesting and general case where the effective potential V{4') considerably changes during inflation. This 
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is the case, for example, in the theories with potentials and e°"^ In such theories one has to treat the boundary 
conditions much more carefully in order to find the correct rate of the volume growth. 

Another approach is to investigate the probability distribution Pp(0,i), which takes into account the exponential 
growth of the volume of domains filled by the inflaton field (f> 0. Solutions for this probability distribution were 
first examined in for the case of chaotic inflation in the theories with potentials 0". It was shown that if 

the initial value of the scalar field is greater than some critical value cj)* , then the probability distribution Pp{(j3,t) 
permanently moves to greater and greater fields (f>. This process continues until the maximum of the distribution 
Pp{(j),t) approaches the field (f>p, at which the effective potential of the field becomes of the order of Planck density 
Mp, where the standard methods of quantum field theory in a curved classical space are no longer valid. By the 
methods used in [^Il0l| it was impossible to check whether Pp{(f>, t) asymptotically approaches any stationary regime 
in the classical domain <j) < <j)p. 

More generally, one may consider the probability distribution which shows the fraction of volume of the 

Universe filled by the field 4> at the time if originally, at the time t = 0, the part of the Universe which we study 
was filled by the field 0o — X- Several important steps towards the investigation of this probability distribution in 
the theories </>" were made by Nambu and Sasaki [23 (see also [2^ E3) ^"^^ by Mijic 0|. Their papers contain 
many beautiful insights, and we will use many results obtained by these authors. However, Mijic [41| did not have 
a purpose in obtaining a complete expression for the stationary distribution Pp{4>, t). The corresponding expressions 
were obtained for various potentials in jsl]. Unfortunately, the stationary distribution P„( 6, t) obtained in 

was almost entirely concentrated at (f)':^ 4>p, i.e. at V{4)) 3> Mp, where the methods used in [23 were inapplicable. 

In the present paper we will argue that the self-reproduction of the inflationary Universe effectively kills itself at 
densities approaching the Planck density. This leads to the existence of a stationary probability distribution Pp{4i, t) 
concentrated entirely at sub-Planckian densities V{4)) < Mp in a wide class of theories leading to chaotic inflation 

The new cosmological paradigm based on inflationary cosmology is very unusual. Instead of the Universe looking 
like an expanding ball of fire, we are considering now a huge fractal consisting of many inflating balls producing new 
balls, producing new balls, etc. In order to make the new cosmological concepts more visually clear, we include in 
this paper a series of figures which show the results of computer simulations of stochastic processes in the inflationary 
Universe |43| . 

The plan of the paper is the following. In Section ^1 we give a short review of chaotic inflation. We discuss the 
generation of density perturbations and explain how they lead to the process of self-reproduction of the Universe. 
In Section imi we discuss the interrelations between stochastic and Euclidean approaches to quantum cosmology and 
explain the main advantages of the probability distribution Pp ((/>, t) for the description of the global structure of 
inflationary Universe. Section llVI brieflv describes the results of the computer simulation of stochastic processes in 
inflationary Universe. Section IVI contains the basics of our approach to finding a stationary probability distribution 
Pp{(j),t\x)- In Section we describe the analogous investigation for a different parametrization of time. Section IVlII 
contains a general discussion of the consequences of our results, which are summarized in a short concluding Section 

EnH 

To simplify our notation, throughout this paper we will use the system of units where Mp = 1. 



II. SELF-REPRODUCING CHAOTIC INFLATIONARY UNIVERSE 



A. Chaotic Inflation 



We will begin our discussion of chaotic inflation with the simplest model based on the theory of a scalar field 
minimally coupled to gravity, with the Lagrangian 

L=-^R+ld,<j>d'^cj,-V{<l)) . (1) 

f DTT 2 

Here G = M^'^ = 1 is the gravitational constant, R is the curvature scalar, and ¥{((>) is the effective potential of the 
scalar field. If the classical field cj) is sufficiently homogeneous in some domain of the Universe (see below), then its 
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behavior inside this domain is governed by the equations 

'<j) + 3i/0 = -dV/d(j) , (2) 

Here H — a/a, a{t) is the scale factor of the Universe, k ~ +1, —1, or for a closed, open or flat Universe, respectively. 

The simplest version of the theory is the theory of a massive noninteracting scalar field with the effective 
potential V{<j)) = ^m'^(j>'^, where m is the mass of the scalar field cj), m ^ 1. If the field (/> initially is sufficiently large, 
0^1, then one can show that the functions (f)(t) and a(t) rapidly approach the asymptotic regime 

^(t)^0o-^^ (4) 

a(t) =aoexp (2^(02 _ 02 (t))) . (5) 

Note that in this regime the second derivative of the scalar field in eq. JSJ and the term fc/a^ in ^ can be neglected. 
The last statement means that the Universe becomes locally flat. 

According to |0J and (0), during the time r ^ (f)/ni the relative change of the field remains small, the effective 
potential V{(f>) changes very slowly and the Universe expands quasi-exponentially: 

a{t + At) - a{t) cxp{HAt) (6) 

for At < T ^ (f>/m. Here 

H = 2^m^. (7) 

Note that r > H^^ for > f . 

The regime of the slow rolling of the field and the quasi-exponential expansion (inflation) of the Universe ends 
at ^ 0e. In the theory under consideration, 0e ^ 0.2. At < 0e the field oscillates rapidly, and if this field 

interacts with other matter fields (which are not written explicitly in eq. JQl), its potential energy V{(/)) ~ ^ ^ 
is transformed into heat. The reheating temperature T/j may be of the order m^^'^ or somewhat smaller, depending 
on the strength of the interaction of the field with other fields. It is important that Tr does not depend on the 
initial value 0o of the field 0. The only parameter which depends on 0o is the scale factor a{t), which grows e'^^'^" 
times during inflation. 

All results obtained above can be easily generalized for the theories with more complicated effective potentials. For 
example, during inflation in the theories with V{4>) = - 0" one has 



For all n. 



Inflation ends at 



'/'^W=0o"^-^y^^ forn^4, (8) 



0(t) =00 exp forn = 4. (9) 



a(i)=aoexp(^(02-02(t))) . (10) 
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Note, that in all realistic models of elementary particles spontaneous symmetry breaking occurs on a scale which 
is many orders of magnitude smaller than 1 (i.e. smaller than Mp). Therefore all results which we obtained remain 
valid in all theories which have potentials V{4>) ^ ^ 0" at 9!) ^ 1, independently of the issue of spontaneous symmetry 
breaking which may occur in such theories at small (j). 

As was first pointed out in |23|, chaotic inflation occurs as well in the theories with exponential effective potentials 
V{4)) = Vb e""^ for sufficiently small a. It was shown later |^ that in the theories with 

V{<P) = e"* (12) 

the Einstein equations and equations for the scalar field have an exact solution: 

(j){t) = ^o-- In^ , (13) 



a(t) = aot^ P= — ■ (14) 

Note that here we are dealing with the power law expansion of the Universe. It can be called inflation if p 3> 1, which 
implies that a <C VIStt ^ 7 [23 . In this case the Hubble constant H almost docs not change within the Hubble time 
~ H"^, and expansion looks quasiexponential, like in cq. ©. Inflation in the theory with the exponential potential 
never ends; it occurs at all 0. In order to make this theory realistic, one should assume that at small (/) the effective 
potential becomes steep and inflation ends a,i (j) < (pe- Without any loss of generality one can assume that (f)^ = 0. In 
this case the parameter Vb gives the value of the effective potential at the end of inflation. 



B. Initial conditions for inflation 



Let us consider first a closed Universe of initial size / ~ 1 (in Planck units), which emerges from the space-time foam, 
or from singularity, or from 'nothing' in a state with the Planck density p ^ \. Only starting from this moment, i.e. at 
p < 1, can we describe this domain as a classical Universe. Thus, at this initial moment the sum of the kinetic energy 

density, gradient energy density, and the potential energy density is of the order unity: ^(jp' + i(i9i0)^ + V{4>) ~ 1. 

We wish to emphasize, that there are no a priori constraints on the initial value of the scalar field in this domain, 
except for the constraint + ^(9^0)'^ + ^(0) ^ 1- Lot us consider for a moment a theory with V{(t)) = const. 
This theory is invariant under the shift (j) ^ (f) + a. Therefore, in such a theory all initial values of the homogeneous 
component of the scalar field (j) are equally probable. Note, that this expectation would be incorrect if the scalar field 
should vanish at the boundaries of the original domain. Then the constraint ^{di(j))^ <^ 1 would tell us that the scalar 
field cannot be greater than 1 inside a domain of initial size 1. However, if the original domain is a closed Universe, 
then it has no boundaries. (We will discuss a more general case shortly.) 

The only constraint on the average amplitude of the field appears if the effective potential is not constant, but 
grows and becomes greater than the Planck density at > (j)p, where V{(f>p) = 1. This constraint implies that (j) <; (j)p, 
but it does not give any reason to expect that (f> <^ (j)p. This suggests that the typical initial value (j>o of the field cf) in 
such a theory is 

(/)o ~ 0p . (15) 

Thus, we expect that typical initial conditions correspond to ^ip^ ^ ^{dicj))^ ^ V{(f>) ~ 0{1). Note that if by 

any chance ^(p^ + ^{dicj))'^ <^ ^{4') in the domain under consideration, then inflation begins, and within the Planck 

time the terms ■^(p^ and -^(9^0)^ become much smaller than V((f>), which ensures continuation of inflation. It seems 
therefore that chaotic inflation occurs under rather natural initial conditions, if it can begin at V{<l>) ~ 1 llllllfl]. 

The assumption that inflation may begin at a very large has important implications. For example, in the theory 
one has 

^o^^p^ . (16) 
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Let us consider for definiteness a closed Universe of a typical initial size 0(1). Then, according to lO, the total size 
of the Universe after inflation becomes equal to 

I ^ exp (27r(/.2) ^ exp . (17) 

For m 10~^ (which is necessary to produce density perturbations ^ ^ 10~^, see below) 

I - exp (27rl0i2) > lQ^^^\m . (18) 

Thus, according to this estimate, the smallest possible domain of the Universe of initial size 0{Mp^) ^ IQ^^^cm after 
inflation becomes much larger than the size of the observable part of the Universe ^ lO^^cm. This is the reason why 
our part of the Universe looks flat, homogeneous and isotropic. The same mechanism solves also the horizon problem. 
Indeed, any domain of the Planck size, which becomes causally connected within the Planck time, gives rise to the 
part of the Universe which is much larger than the part which we can see now. 

In what follows we will return many times to our conclusion that the most probable initial value of the scalar field 
corresponds to ~ 4>p. There were many objections to it. Even though all these objections were answered many 
years ago |7[ |23| , we need to discuss here one of these objections again. This is important for a proper understanding 
of the new picture of the evolution of the Universe in inflationary cosmology. 

Assume that the Universe is not closed but infinite, or at least extremely large from the very beginning. (This 
objection does not apply to the closed Universe scenario discussed above.) In this case one could argue that our 
expectation that (/)o ~ ^ 1 is not very natural ji^. Indeed, the conditions (dicj))'^ <^ 1 and (po ^ (j)p ^ I imply 
that the field </> should be of the same order of magnitude ~ (/ip 1 on a length scale at least as large as (ftp, which 
is much larger than the scale of horizon Z ~ 1 at the Planck time. But this is highly improbable, since initially (i.e., 
at the Planck time ) there should be no correlation between values of the field (p in different regions of the Universe 
separated from one another by distances greater than 1. The existence of such correlation would violate causality. As 
it is written in f43| , the scalar field (j) must be smooth on a scale much greater than the scale of the horizon, which 
does not sound very chaotic. 

The answer to this objection is very simple [SIlOl- We have absolutely no reason to expect that the overall energy 
density p simultaneously becomes smaller than the Planck energy density in all causally disconnected regions of an 
infinite Universe, since that would imply the existence of an acausal correlation between values of p in different 
domains of Planckian size Ip I. Thus, each such domain at the Planck time after its creation looks like an isolated 
island of classical space-time, which emerges from the space-time foam independently of other such islands. During 
inflation, each of these islands independently acquires a size many orders of magnitude larger than the size of the 
observable part of the Universe. A typical initial size of a domain of classical space-time with p <J 1 is of the order 
of the Planck length. Outside each of these domains the condition p <^ 1 no longer holds, and there is no correlation 
between values of the field 4> in different disconnected regions of classical space-time of size 1. But such correlation 
is not necessary at all for the realization of the inflationary Universe scenario: according to the 'no hair' theorem for 
de Sitter space, a sufficient condition for the existence of an inflationary region of the Universe is that inflation takes 
place inside a region whose size is of order H^^. In our case this condition is satisfied. 

We wish to emphasize once again that the confusion discussed above, involving the correlation between values of 
the field 4> in different causally disconnected regions of the Universe, is rooted in the familiar notion of a very large 
Universe that is instantaneously created from a singular state with p = oo, and instantaneously passes through a state 
with the Planck density p — I. The lack of justification for such a notion is the very essence of the horizon problem. 
Now, having disposed of the horizon problem with the aid of the inflationary Universe scenario, we can possibly 
manage to familiarize ourselves with a different picture of the Universe. In this picture the simultaneous creation 
of the whole Universe is possible only if its initial size is of the order 1, in which case no long-range correlations 
appear. Initial conditions should be formulated at the Planck time and on the Planck scale. Within each Planck-size 
island of the classical space-time, the initial spatial distribution of the scalar field cannot be very irregular due to the 
constraint (dicj))'^ < 1. But this does not impose any constraints on the average values of the scalar field (j) in each of 
such domains. One should examine all possible values of the field (f) and check whether they could lead to inflation. 

Note, that the arguments given above ^ suggest that initial conditions for inflation are quite natural only if 
inflation begins as close as possible to the Planck density. These arguments do not give any support to the models 
where inflation is possible only at densities much smaller than 1. And indeed, an investigation of this question shows, 
for example, that a typical closed Universe where inflation is possible only at U((/)) ^ 1 collapses before inflation 
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begins. Thus, inflationary models of that type require fine-tuned initial conditions [46| . and apparently cannot solve 
the flatness problem. 

Does this mean that we should forget all models where inflation may occur only at V{(j)) ^1? Does this mean, in 
particular, that the 'natural inflation' model is, in fact, absolutely unnatural? 

As we will argue in the last section of this paper, it may be possible to rescue such models either by including 
them into a more general cosmological scenario ('hybrid inflation' |27|'). or by using some ideas of quantum cosmology, 
which we are going to elaborate. 



C. Quantum Fluctuations and Density Perturbations 



According to quantum fleld theory, empty space is not entirely empty. It is filled with quantum fiuctuations of all 
types of physical fields. These fluctuations can be considered as waves of physical flelds with all possible wavelengths, 
moving in all possible directions. If the values of these fields, averaged over some macroscopically large time, vanish, 
then the space filled with these fields seems to us empty and can be called the vacuum. 

In the exponentially expanding Universe the vacuum structure is much more complicated. The wavelengths of all 
vacuum fluctuations of the scalar field (p grow exponentially in the expanding Universe. When the wavelength of any 
particular fluctuation becomes greater than H~^, this fluctuation stops oscillating, and its amplitude freezes at some 
nonzero value 6(j){x) because of the large friction term SHcj) in the equation of motion of the field 0. The amplitude 
of this fluctuation then remains almost unchanged for a very long time, whereas its wavelength grows exponentially. 
Therefore, the appearance of such a frozen fluctuation is equivalent to the appearance of a classical field 5(j){x) that 
does not vanish after averaging over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the creation of more and more new 
perturbations of the classical field with wavelengths greater than . The average amplitude of such perturbations 
generated during a time interval H^^ (in which the Universe expands by a factor of e) is given by^ 

mx)\ « ^ . (19) 

If the field is massless, the amplitude of each frozen wave does not change in time at all. On the other hand, phases of 
each waves are random. Therefore, the sum of all waves at a given point fluctuates and experiences Brownian jumps 
in all directions. As a result, the values of the scalar field in different points become different from each other, and 
the corresponding variance grows as 

{^') - £ ^ (20) 

i.e. 

H 



VW) = ^VHt. (21) 

This result, which was first obtained in "J?], is an obvious consequence of eq. (|19|l . if one considers the Brownian 
motion of the field cj). One should just remember that the field makes a step ^ each time At = H^^, and the total 
number of such steps during the time t is given by iV = Ht. 

If the field (p is massive, with m <^ H , then the amplitudes of the fiuctuations frozen at each time interval At ~ 
are given by the same equation as for the massless field, but later the amplitudes of long- wavelength fiuctuations begin 
to decrease slowly and the variance of the long- wavelength fiuctuations, instead of indefinite growth H21|l . approaches 
the limit ^ 



^ To be more precise, the amplitude of each wave is greater by \/2, but this factor disappears after taking an average over all phases. 
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If perturbations grow while the mean value of the field decreases, the behavior of the variance becomes a little bit 
more complicated, see Section Hi Ul 

Fluctuations of the field (f> lead to adiabatic density perturbations 

Sp ^ V'{^)S(t> , (23) 

which grow after inflation, and at the stage of the cold matter dominance acquire the amplitude 

Sp 48/2^1^ 

p " 5 V 3 V'icj)) ■ ^ > 

Here is the value of the classical field (t){t) (4), at which the fluctuation we consider has the wavelength I ^ ^ 

rn_ 
2 



H ^{(j)) and becomes frozen in amplitude. In the theory of the massive scalar field with V{(j)) — ^-^'^-'^ 



Sp 24 /tt _ 2 



-mr . (25) 

Taking into account Q), jSJ and also the expansion of the Universe by about 10'^'^ times after the end of inflation, 
one can obtain the following result for the density perturbations with the wavelength l{cm) at the moment when these 
perturbations begin growing and the process of the galaxy formation starts: 

— - 0.8 m In Z(cm) . (26) 
P 

This implies that the density perturbations acquire the necessary amplitude y ~ 10^^ on the galaxy scale, Ig ~ 10^^ 
cm, if m ~ 10~^, in Planck units, which is equivalent to 10^'^ GeV. 

Similar constraints on the parameters of inflationary models can be obtained in all other theories discussed in 
the previous subsection. We will mention here the constraints on the theory with exponential effective potential 
¥{(/)) — Vq e°"^. Density perturbations in this theory have the following dependence on a, Vq and the length scale I 
measured in cm: 

--15 I ; p= — (27) 

pa 

This spectrum in realistic models of inflation should not significantly (much more than by a factor of 2) change 
between the galaxy scale Ig ^ 10^^ cm and the scale of horizon Ih ^ 10^* cm, and should be of the order of 10^^ on 
the horizon scale. This leads to the constraints a <^ I and Vq ^ lO^^^a^. 

In our previous investigation we only considered local properties of the inflationary Universe, which was quite 
sufficient for description of the observable part of the Universe of the present size I ~ 10^^ cm. For example, in the 
model of a massive field <p with m ^ 10~^ our Universe, in accordance with (|26l) . remains relatively homogeneous up 
to the scale 

l<l* e^/"'cm - 10^-^°' cm . (28) 

The density perturbations on the scale /* were formed at the time when the scalar field (j){t) was equal to 0*, where 
(see ^) 

^ 5 X 10^ , (29) 



Note that V{(j)*) ~ ^(0*)^ ~ m <C 1. On scales I > I* the Universe becomes extremely inhomogeneous due to 
quantum fluctuations produced during inflation. 

Similar conclusion proves to be correct for all other models discussed above. In the models with V{(j)) = -^0" the 
critical value of the field 6 can be estimated as 
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whereas in the models with the exponential potential 



The corresponding scale I* — l{4)*) gives us the typical scale at which the Universe still remains relatively homogeneous 
and can be described as a Friedmann Universe. In all realistic models this scale is much larger than the size of the 
observable part of the Universe, but much smaller than the naive classical estimate of the size of a homogeneous part 
of the Universe H17|l . 

We are coming to a paradoxical conclusion that the global properties of the inflationary Universe are determined 
not by classical but by quantum effects. Let us try to understand the origin of such a behavior of the inflationary 
Universe. 



D. Self- Reproducing Universe 



A very unusual feature of the inflationary Universe is that the processes separated by distances I greater than H^^ 
proceed independently of one another. This is so because during exponential expansion the distance between any two 
objects separated by more than is growing with a speed v exceeding the speed of light. As a result, an observer 
in the inflationary Universe can see only the processes occurring inside the horizon of the radius H^^. 

An important consequence of this general result is that the process of inflation in any spatial domain of radius H^^ 
occurs independently of any events outside it. In this sense any inflationary domain of initial radius exceeding 
can be considered as a separate mini- Universe, expanding independently of what occurs outside it. This is the essence 
of the "no-hair" theorem for de Sitter space, which we already mentioned in Section 2.2. 

To investigate the behavior of such a mini-Universe, with an account taken of quantum fluctuations, let us consider 
an inflationary domain of initial radius H^^ containing sufficiently homogeneous field with initial value 0^1 (assume 

V{4)) = ^ for simphcity). Equation Q tells us that during a typical time interval At = iJ^^ the field inside this 
domain will be reduced by 

By comparison of (|19|l and l|32|l one can easily see that if is much less than </>* ~ , then the decrease of the field 
due to its classical motion is much greater than the average amplitude of the quantum fluctuations generated 
during the same time. But for greater (up to the classical limit of about 10^), (50(a:) will exceed A(/), i.e. the 
Brownian motion of the field (/) will become more rapid than its classical motion. Because the typical wavelength 
of the fluctuations i50(a;) generated during this time is ^ the whole domain after At effectively becomes divided 
into separate domains (mini-Universes) of radius each containing almost homogeneous field (/) — A(/) -I- <50. 

We will call these domains "/i-regions" [s^, 113 , to indicate that each of them has the radius coinciding with the 
radius of the event horizon . In almost half of these domains (i.e. in e'^/2 ~ 10 /i-regions) the field ^ grows by 
|(50(x)| — A(/) « |(5(/)(a:)| — i?/27r, rather than decreases. During the next time interval At = H^^ the field grows again 
in the half of the new /i-regions. Thus, the total number of /i-regions containing growing field becomes equal to 
(e"^/2)^ = e^^^"'"^^). This means that until the fluctuations of fleld grow sufficiently large, the total physical volume 
occupied by permanently growing field (i.e. the total number of /i-regions containing the growing field 0) increases 
with time like exp[(3 — In 2) i7t]. The growth of the volume of the Universe at later stages becomes even faster. 

For example, let us consider those /i-regions where the field grows permanently, i.e. where the jumps of the field </) 
are always positive, (50 ~ H/2-k. Since this process occurs each time H"^ , the average speed of growth of the scalar 
field in such domains is given by 

d4> _ H^ct>) _ 4U(0) 
dt " 27r " 3 ■ 

In the theory of a massive noninteracting field with V{(f)) = m^(/)^/2 the solution of this equation is 

rHO^^o-^-^, (34) 
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where (j)Q is the initial value of the field. Thus, within the time 



t = r. n , (35) 

the field (f> in those domains becomes infinitely large. 

Of course, this solution cannot be trusted when the field approaches the value (j)p ~ ^ 1, where the effective 
potential becomes of the order of one (i.e. when the energy density becomes comparable with the Planck density). 
In fact, we are going to argue in Section FV Al that the growth of the field (p stops when it approaches <j)p- What is 
important though is that within a finite time (|35|l a finite part of the volume of the Universe approaches a state of 
the maximally high energy density. After that moment, a finite portion of the total volume of the Universe will stay 
in a state with the Planck density and will expand with the maximal possible speed, a{t) ~ e*^*, where C = 0(1) 
(i.e. C ^ Mp). All other parts of the Universe will expand much slower. Consequently, the parts of the Universe with 
(j) ^ (f)p will give the main contribution to the growth of the total volume of the Universe, and this volume will grow 
as fast as if the whole Universe were almost at the Planck density. 

Thus, if these considerations are correct, the main part of the physical volume of the Universe is the result of the 
expansion of domains with nearly the maximal possible field value, (j) ^ 4'p, for which V{4>) is close to the Planck 
boundary V ^ 1. There are also exponentially many domains with smaller values of 4>. Those domains, in which (jj 
eventually becomes sufficiently small, give rise to the mini-Universes of our type. In such domains, (j) eventually rolls 
down to the minimum of V{(j)), and these mini-Universes are subsequently describable by the usual Big Bang theory. 
However a considerable part of the physical volume of the entire Universe remains forever in the inflationary phase 

i 

Let us return to the discussion of the properties of our stochastic processes. Eq. (|35|l suggests that it takes time 
t ~ 2rW-^4>a ^^til the Planck density domains will dominate the speed of the growth of the volume of the Universe. To 
be sure that this conclusion is correct, one should take into account that at the beginning of the Planckian expansion 
of the domains where the field (j) was permanently growing, their initial volume could be somewhat smaller than the 
total volume of all other domains, where the field (f> was not permanently growing. Indeed, only in a small portion of 
the original space the scalar field jumps up permanently; typically it will jump in both directions. However, one can 
show that for m ^ 1 the corresponding corrections are small. One should check also that this result is not modified 
by the possibility of larger quantum jumps of the scalar field. Indeed, the typical amplitude of the quantum jump 
within the time is given hy 5(j} ^ However, there exists a nonvanishing probability that the field within 
the same time will jump much higher. Even though this probability is exponentially suppressed, such domains will 
start their expansion with the Planckian speed earlier, and their relative contribution to the volume of the Universe 
could become significant. 

Assuming that this probability of a jump — > -f (50 is given by the Gaussian distribution with variance one 
obtains 

P(.0)=exp(-^). (36) 



For example, the probability of a jump from (po 4> to ip can be estimated by 

bo) 

For the theory of a massive scalar field this gives the following probability of a jump to 4> = 4>p: 

4>l 

Even though the volume of this part of the Universe grows with the Planckian speed, a ~ e^^*, c ~ 0(1), it takes 
longer than t ~ J^'" ^ for this part of the Universe to grow to the same size as the rest of the Universe. This time is 

much larger than the time t = ^^2^^ which we obtained neglecting the possibility of large jumps. 

In fact, our estimate of the probability of large jumps was even too optimistic. Indeed, if the jump is very high, 
then the gradient energy of the corresponding quantum fluctuation on a scale is given by jj^^^- This quantity 
may become greater than the potential energy density of the scalar field V {(()). In such a case the domain where the 
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jump occurs does not inflate. To overcome this difficulty one should consider jumps on a larger length scale, which 
further reduces their probability to 

2 J 4 

P{4>) ^ exp(-c^) , (39) 

where c = 0(1) |49| |. One could conclude, therefore, that large quantum jumps are really unimportant, and the main 
features of the process can be understood neglecting the possibility of jumps greater than ^ . 

However, one should keep in mind that eq. is valid only for 3> 0*. For < 0*, the mean decrease 

of the field (p during the typical time H^^ is always more significant than the typical amplitude of the jump 
Therefore, any process of self-reproduction of the Universe with a small initial value of the field (j) is possible only 
with an account taken of large quantum jumps. The probability of such jumps is exponentially suppressed. Indeed, 
by comparing eqs. (|37l) and 117(1 one can see that the strong suppression of large jumps cannot be compensated even 
by the exponentially large inflation after such a jump. Therefore in a typical /i-region with (j) < (f>* the process of the 
Universe self-reproduction does not occur and inflation eventually ceases to exist. However, importance of large jumps 
should not be overlooked: one such jump may be enough to start an infinite process of the Universe self-reproduction. 
This process may lead to spontaneous formation of an inflationary Universe from the ordinary Minkowski vacuum 
and to a possibility of the Universe creation in a laboratory |4^. We will return to this question in Section 7.2. 

Note that the process of the self-reproduction of the Universe may occur not only in the chaotic inflation scenario, 
but in the old jl^ and the new inflationary scenarios as well. In the corresponding models some part of the 
volume of the Universe can always remain in a state corresponding to a local extremum of V{(j)) a.t cf) — 0. In the 
chaotic inflation scenario in the theories 0" or e"'^, the results are more surprising. Not only can the Universe stay 
permanently on the top of a hill as in the old and new inflationary scenario; it can also climb perpetually up the wall 
towards the largest possible values of its potential energy density |^. 



III. STOCHASTIC APPROACH TO INFLATION AND THE WAVE FUNCTION OF THE UNIVERSE 



A. Euclidean Approach 



Now we would like to compare our methods with other approaches to quantum cosmology. One of the most ambitious 
approaches to cosmology is based on the investigation of the Wheeler-DeWitt equation for the wave function "if of 
the Universe [sj. However, this equation has many different solutions, and a priori it is not quite clear which one of 
these solutions describes our Universe. 

A very interesting idea was suggested by Hartle and Hawking [s^- According to their work, the wave function of 
the ground state of the Universe with a scale factor a filled with a scalar field cj) in the semi-classical approximation 
is given by 

*o(a,0)-exp(-5£(a,0)) . (40) 

Here S'£;(a, 0) is the Euclidean action corresponding to the Euclidean solutions of the Lagrange equation for air) 
and ^(r) with the boundary conditions a(0) — a, (l>(0) = 4>. The reason for choosing this particular solution of the 
Wheeler-DeWitt equation was explained as follows [33|. Let us consider the Green's function of a particle which 
moves from the point (0, t') to the point x, t: 

<x,0|0,t' >=^^'„(x)«'„(0)exp {iEnit^t' j) 

71 

= J dx{t)exp{iS{x{t))) , (41) 

where "^fn is a complete set of energy eigenstates corresponding to the energies E„ > 0. 

To obtain an expression for the ground-state wave function ^o(x), one should make a rotation t —it and take 
the limit as r ^ — oo. In the summation (|41|l only the term n = with the lowest eigenvalue i?o — survives, and 
the integral transforms into J (ix(r) exp(— S'£;(r)). Hartle and Hawking have argued that the generalization of this 
result to the case of interest in the semiclassical approximation would yield (|40|l . 
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The gravitational action corresponding to the EucHdean section S4 of de Sitter space dSi with a{T) = H ^ (cf)) cos Ht 
is negative, 



SE{a,^) = -lJdv[{^y-a' + -^ 



Svr 
T 



mv{(i)) 



(42) 



Here rj is the conformal time, rj = J ^4|j, A = 8nV. Therefore, according to 



*o(a, 4>) exp (^-Ssia, 4>) 



exp 



V 16^(0) 



(43) 



This means that the probabiHty P of finding the Universe in the state with cf) = const, a = H ^{4>) = y^^^j^ is 
given by 

P(^)H^oP-exp(^) . (44) 



This expression has a very sharp maximum as V{(l)) — > 0. Therefore the probability of finding the Universe in a 
state with a large field (j) and having a long stage of inflation becomes strongly diminished. Some authors consider 
this as a strong argument against the Hartle-Hawking wave function. However, nothing in the 'derivation' of this 
wave function tells that it describes initial conditions for inflation; in fact, in the only case where it was possible to 
obtain eq. 143() by an alternative method, the interpretation of this result was quite different, see Section 3.2. 

There exists an alternative choice of the wave function of the Universe. It can be argued that the analogy between 
the standard theory (|41|l and the gravitational theory (|42|l is incomplete. Indeed, there is an overall minus sign in 
the expression for SEia,(j)) H42|) . which indicates that the gravitational energy associated with the scale factor a is 
negative. (This is related to the well-known fact that the total energy of a closed Universe is zero, being a sum of 
the positive energy of matter and the negative energy of the scale factor a.) In such a case, to suppress terms with 
En < Eq and to obtain from (|42|) one should rotate t not to —it, but to +iT. This gives ^3 

^-0(0, ~ exp(-|5i5(a, 4>)\) ~ exp j , (45) 

and 

P(</.)~|*o(a,0)P~exp(-2|5i,(a,0)|) ^exp^-^^] . (46) 



An obvious objection against this result is that it may be incorrect to use different ways of rotating t for quantization 
of the scale factor and of the scalar field. However, the idea that quantization of matter coupled t o g ravity can be 
accomplished just by a proper choice of a complex contour of integration, though very appealing |50|, may be too 
optimistic. We know, for example, that despite many attempts to suggest Euclidean formulation (or just any simple 
set of Feynman rules) for nonequilibrium quantum statistics or for the field theory in a nonstationary background, 
such formulation does not exist yet. It is quite clear from H41|l that the t —it trick would not work if the spectrum 
En were not bounded from below. Absence of equilibrium, of any simple stationary ground state seems to be a typical 
situation in quantum cosmology. In some cases where a stationary or quasistationary ground state does exist, eq. H44|l 
may be correct, see Section 3.2. In a more general situation it may be very difficult to obtain any simple expression for 
the wave function of the Universe. However, in certain limiting cases this problem is relatively simple. For example, 
at present the scale factor a is very big and it changes very slowly, so one can consider it to be a C-number and 
quantize matter ffelds only by rotating t — s- — ir. On the other hand, in the inffationary Universe the evolution of the 
scalar field is very slow; during the typical time intervals 0{H~^) it behaves essentially as a classical field, so one can 
describe the process of the creation of an inffationary Universe filled with a homogeneous scalar field by quantizing 
the scale factor a only and by rotating t — + it. Eq. which was first obtained in by the method described 

above, later was obtained also by another method by Zeldovich and Starobinsky Rubakov j^^, and Vilenkin [T^ . 
This result can be interpreted as the probability of quantum tunneling of the Universe from a = (from "nothing" ) 
to a = H~^{(j>). Therefore the wave function H45|l is called 'tunneling wave function'. In complete agreement with our 
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previous argument, eq. (|46() predicts that a typical initial value of the field (j) is given by V{(j)) ~ 1 (if one does not 
speculate about the possibility that V{(l)) ^ 1), which leads to a very long stage of inflation. 

Unfortunately, there is no rigorous derivation of either l|44(l or H46|l . and the physical meaning of creation of 
everything from 'nothing' is far from being clear. However, one may give a simple physical argument explaining that 
(|46|l has a better chance to describe the probability of quantum creation of the Universe than (|44|l [T^ . 

Examine a closed de Sitter space with energy density V{(j)). Its size behaves as H^^ cosh Ht. Thus, its minimal 
volume at the epoch of maximum contraction is of order ~ V~'^^'^, and the total energy of the scalar field 

contained in de Sitter space at that instant is approximately V x V^"^^^ ~ V^^^^. Thus, to create the Universe with 
the Planckian energy density V{4>) ~ 1 one needs only a Planckian fluctuation of energy AE ^ 1 on the Planck scale 

^ 1, whereas to create the Universe with V{(j)) ^ 1 one needs a very large fluctuation of energy AE ~ 
on a large scale R-^ - y~i/2. This means that, in accordance with (|46|l . the probability of quantum creation of the 
Universe with V{(f)) <C 1 should be strongly suppressed, whereas no such suppression is expected for the probability 
of creation of an inflationary Universe with V{(j)) ^ 1. 

This argument suggests that, in agreement with our discussion of initial conditions in Section 2.2, inflation appears 
in a much more natural way in the theories where it is possible for V{4>) ~ 1 (e.g. in the theories with the effective 
potentials 0" or e"*^). rather than in the theories where it may occur only at V{4>) <C 1 (theories of this class include 
'natural inflation' [23, 'hyperextended inflation' 51], etc.). 

A deeper understanding of the physical processes in the inflationary Universe is necessary in order to investigate the 
wave function of the Universe \l/o(a, 0), to suggest a correct interpretation of this wave function and to understand 
a possible relation of this wave function to our results concerning the process of self-reproduction of inflationary 
Universe. With this purpose we will try to investigate the global structure of the inflationary Universe, and go beyond 
the minisuperspace a ppr oach used in the derivation of (|40|) and (|46|l . This can be done with the help of a stochastic 
approach to inflation |lCtl34j. which is a more formal way to investigate the Brownian motion of the scalar field (j). As 
we will see, this approach may somewhat modify our conclusions concerning the possibility to realize various models 
of inflation. 



The evolution of the fluctuating field (f) in any given domain can be described with the help of its distribution 
function P{4>)^ or in terms of its average value (f) in this domain and its variance A = {StjP'). However, one will 
obtain different results depending on the method of averaging. One can consider the probability distribution Pc{(t), t) 
over the non-growing coordinate volume of the domain (i.e. over its physical volume at some initial moment of 
inflation). This is equivalent to the probability to find a given field at a given time at a given point.^ Alternatively, 
one may consider the distribution Pp{(t>, t) over its physical (proper) volume, which grows exponentially at a different 
rate in different parts of the domain. In many interesting cases the variance Ac of the field cf) in the coordinate volume 
remains much smaller than (pc- In such cases the evolution of the averaged field 0c, can be described approximately 
by ©-©■ However, if one wishes to know the resulting spacetime structure and the distribution of the field (f) after 
(or during) inflation, it is more appropriate to take an average (j)p over the physical volume, and in some cases the 
behavior of (j)p and Ap differs considerably from the behavior of 0c and Ac. 

We will begin with the description of the distribution Pc{4>, t). This is a relatively simple problem. As we learned 
already, the field at a given point behaves like a Brownian particle. The standard way to describe the Brownian 
motion is to use the diffusion equation, where instead of the position of the Brownian particle x one should write 0: 



Here V is the coefficient of diffusion, k is the mobility coefhcient, — ^ is the analog of an external force F. The 
parameter [3 reflects some ambiguity which appears when one derives this equation in the theories with T) depending 



B. Stochastic Approach to Inflation 




(47) 



^ Note, that by the field <j) we understand here its classical long- wavelength component, with the wavelength I > H ^. Thus, this field 
remains (almost) constant over each /i-region. 
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on (j>. The value f3 = 1 corresponds to the Ito version of stochastic approach, whereas /? = 1/2 corresponds to the 
version suggested by Stratonovich. 

The standard way of deriving this equation is based on the Langevin equation, which in our case can be written in 
the following way [3 : 

Here ^(t) is some function representing the effective white noise generated by quantum fluctuations, which leads to the 
Brownian motion of the classical field (f>. The derivation of this equation, finding the function S,(t) and the subsequent 
derivation of the diffusion equation is rather complicated [s^. I3It| . Meanwhile, the main purpose of the derivation is 
to establish the functional form of V and n. This can be done in a very simple and intuitive way 

Indeed, let us study two limiting cases. The first case is a classical motion without any quantum corrections. In 

this case the field (j> during inflation satisfies equation (|48|l without the last term; see eq. (01 and discussion after it. 

Comparing this equation and the standard definition of the mobility coefficient (i = kF) one easily establishes that 
1 

mm- 



m our case k = ^ 



Determination of the diffusion coefficient is slightly more complicated but still quite elementary. This coefficient 
describes the process of adding (and freezing) new long-wavelength perturbations. As we mentioned already, the speed 
of this process does not depend on m at m? <C -ff^. Therefore we may get a correct expression for D by investigating 
the theory with a flat effective potential, where the last term in eq. l|T7|l disappears. According to eq. ((^OJ), in such a 
theory 

d 

i - h ■ <^») 

On the other hand, eq. (|47|l without the last term gives 

= 2V I d(P Pc = 2V . 



When going from the first to the second line of this equation, we used the fact that V — const in the theory with the 
flat potential. From equations ^ and lO it follows that V = ^ = '^^J^^ ■ This gives us the diffusion equation 

or, in an expanded form, 

^ ^ ^ (•'*-""^<*) (^'-(*)-=<*.')) - iai-^<^^'.) ■ (-) 

This equation for the case -ff (</>) = const was first derived by Starobinsky 1341 : for a more detailed derivation see 
[ssj . ^6]. For the special case = this equation was obtained by Vilenkin |3(]|. This equation can be represented 
in the following useful form, which reveals its physical meaning: 

dPc dJc , , 

where the probability current is given by 



(54) 



Equation (|53|l can be interpreted as a continuity equation, which follows from the conservation of probability. We 
will use this representation in our discussion of the boundary conditions for the diffusion equations. 
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We still did not specify the value of the parameter (3 in this equation. The uncertainty in /? is related to the 

^3/2 / - ■, 

precise definition of the white noise term — ' ^{t) in the Langevin equation. At H — const (de Sitter space) there 
is no ambiguity in definition of the white noise S,{t). However, if one wishes to take into account the dependence 
of H on 4>(t), one may include some part of this dependence into the definition of white noise. Alternatively (and 
equivalently) , one may use several different ways to obtain the diffusion equation from the Langevin equation, which 
also leads to the same uncertainty. (A similar ambiguity is known to appear at the level of operator ordering in the 
Wheeler-DeWitt equation.) In what follows we will use the diffusion equation in the Stratonovich form (/? = 1/2), 

5Pe(0,f) d / ,3/4,.^ 5 ^.3/4,.^„,. , 3r(^) 



dt 3V3^ 



but one should keep in mind that most of the results which we will obtain are not very sensitive to a particular choice 
of (3. For example, a simplest stationary solution [dPc/dt = 0) of equation H52(l would be [23| 



P,(^)^F-3/'/2(^).exp(^- 



The whole dependence on (3 here is concentrated in the subexponential factor V~''^^^'^{4>). 



(56) 



In fact, one can make another step and derive an equation describing the probability distribution P{(j), t\x) that 
the value of the field (t){t) initially (i.e. at t = 0) was equal to x- The derivation of a similar equation (which is called 
'backward Kolmogorov equation') will be contained in Section^ Now we are just presenting the result: 

Note that in this equation one considers </> as a constant, and finds the time dependence of the probability that this 
value of the scalar field was produced by the diffusion from its initial value x during the time t. The simplest 
stationary solution of both eq. (|55|) and eq. (|57|l (subexponential factors being omitted) would be 

^^(^'*l^)-^^p(84))-^^p(-84)) ■ 



This function is extremely interesting. Indeed, the first term in H58|l is equal to the square of the Hartle-Hawking 
wave function of the Universe 14411 . whereas the second one gives the square of the tunneling wave function H46|l ! 

At first glance, this result gives a direct confirmation and a simple physical interpretation of both the Hartle- 
Hawking wave function of the Universe and the tunneling wave function. However, in all realistic cosmological 
theories, in which U(0) = at its minimum, the distributions H56I) . (|58|) are not normalizable. The source of this 
difhculty can be easily understood: any stationary distribution may exist only due to a compensation of a classical 
flow of the field (j) downwards to the minimum of U(0) by the diffusion motion upwards. However, diffusion of the 
field 4) discussed above exists only during inflation, i.e. only for > 1, V{(t)) > V{1) ^ ~ 10~^^ for m ~ 10~^. 
Therefore ()58|l would correctly describe the stationary distribution Pc{4),t\x) in the inflationary Universe only if 
y{4>) ^ 10"^'^ ^ 10*" g-cm^'^ in the absolute minimum of V{(j))^ which is, of course, absolutely unrealistic [Tol |. 

Of course, H58|l is not the most general stationary solution for P^- For example, eq. H52I) in the case /3 = 0, 
V{4)) = ni^2- has a simple constant solution 

Pe(</>,t|x)=-Jc^ . (59) 
m 

Similar solutions were described in [36| : for a very recent discussion see also "5^. However, an interpretation of such 
solutions is obscure, see e.g. 0,0. Indeed, eq. H59(l is a solution of eq. H52|l even in the absence of any diffusion; 
one may simply ignore the first term in eq. 1)52(1 . It may exist only if from the very beginning there were infinitely 
many /i-regions with all possible values of the field (j) equally distributed everywhere from to 00. In this case the 
motion of the whole distribution H59() towards small </) due to the classical rolling of the field does not change the 
distribution. But this presumes that the Universe from the very beginning was infinitely large, that our probability 
distribution was fine-tuned, and that our diffusion equations are valid at infinitely large energy density. There was 
an attempt to interpret these solutions as representing creation of new domains of the Universe from nothing. Such 
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an interpretation could be possible for the distribution Pp, which takes into account increase of the volume of the 
Universe. Meanwhile, the derivation of the diffusion equation for presumes that there is no creation of new space 
in comoving coordinates. For all these reasons (see also below) we will not discuss 'stationary solutions' of the type 
of l|59|) in our paper. 

To get a better understanding of the situation, we will study a general non-stationary solution of H55|) describing 
the time-evolution of the initial distribution Pc{(f>,t = 0) = d{4> — 0o)- This is a rather general form of the initial 
distribution Pf.(4>,t) in a domain of initial size of the order of H^^. Indeed, the typical deviation from homogeneity 
in such domains is given by the amplitude of quantum perturbations on this scale, ^ One can easily check 
that for m ^ 1 this amplitude is always much smaller than (j) during inflation. This explains why the delta-functional 
initial conditions are relevant to our problem. 

The solution of the diffusion equation for the theory with these initial conditions is [To|: 



where (pit) is the slow-rollover solution of the classical equations (0), This equation shows that at the first stage 

of the process, during the time At < "^"^^^ (see eq. (@J), the maximum of the distribution Pc{4>,t) almost does 
not move, whereas the variance grows linearly. Then, at i 3> At, the maximum of Pc{4>,t) moves to = 0, just as 
the classical field 0(t) This shows that in the realistic situations with reasonable initial conditions there are no 
nontrivial stationary solutions of the diffusion equation for Pc'. the field just moves towards the absolute minimum 
of the effective potential and stays there 9]. 

For completeness we shall mention here another solution of H55() . If the initial value of the field (p is very large, 
(f>o > m~^, i.e. if one starts with the spacetime foam with V{(f>o) > 1, then the evolution of the field in the first 
stage (rapid diffusion) becomes more complicated (the naively estimated variance A^ ~ H^t soon becomes greater 
than (/iq). In this case the distribution of the field (p is not Gaussian. The solution of eq. H55|l at the stage of diffusion 
from <^o to some field (j) with V{(f>) ^ 1 is given by 

This solution describes quantum creation of domains of a size / > H^^{4>), which occurs due to the diffusion of the 
field 4> from (j>o ^ m'^ to (j> cpQ. Direct diffusion with formation of a domain filled with the field (p is possible only 

during the time 2c\/37r(/i/m, c = 0(1). At larger times a more rapid process is diffusion to some field (j) > (f) and a 
subsequent classical rolling down from (j) to (p. Therefore one may interpret a distribution Pc{<j)) formed after a time 
At = {2c^/3TT(j)/m) as the probability of quantum creation of a mini-Universe filled with a field (j) |lO|. 

3c \ / 3c 

J ^ ''''P V 4U("0) 

which is in agreement with the previous estimate for the probability of quantum creation of the Universe, eq. 146|) 
|l3 | -|l6 | . However, this result is not a true justification of our expression for the probability of the quantum creation 
of the Universe. First of all, we did not actually determine the constant c. Moreover, the careful analysis shows that 
for the effective potentials steeper than X(p'^ the simple expression H62|l should be modified. For example, for V ^ cj)", 
n > 4, an improved result is [3| 

Now let us try to write an equation of the type of H55() for Pp. Let us first forget about the normalization of the 
distribution Pp. Then the distribution Pp{(j),t) has a meaning of a total number of "points" with a given (j), provided 
that the number of points not only changes due to diffusion, but also grows proportionally to the increase of volume. 
This means that during a small time interval dt the total number of points with the field (f) additionally increases by 
the factor 2>H{(t)) dt. (For (f) — const this would lead to an exponential growth of the number of points corresponding 
to the expansion of the volume ^ e'^^*). This leads to the following equation for the (unnormalized) distribution Pp: 



PA.4>) ~ oxp ( ) ~ exp ( ) , (62) 
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or, in an expanded form, 



This simple derivation was proposed by Zeldovich and one of the present authors immediately after the discovery 
of the self-reproduction of the Universe in the chaotic inflation scenario An alternative (and more rigorous) 

derivation was given in mHa. In Section EH of this paper we will give another derivation of this equation and its 
generalizations, using methods of the theory of branching diffusion processes. 

In the first papers on the self-reproduction of the Universe this equation was not used. Instead of that, it was found 
useful to study solutions of equation II55II and then make some simple estimates which gave a qualitatively correct 
description of the behavior of Pp 0, For example, one can use eq. and take into account that during the 

first period of time At < when the field almost does not move, the volume of domains filled by the field (j) 

increases approximately by exp(3ff ((/))At) ~ exp(ci(/>o(/') where ci = 0(1). After this time the distribution Pp{(t>, At) 
is (approximately) given by 

Ppicj},At) ~ Pc{^,At) ■ exp(3iI(0)At) ^ exp (- ^^^'^4 + ck^o^) , (66) 

where C2 = 0(1). One can easily verify that, if 0o ^ 0* ~ then the maximum of Pp{(j),t) during the time 

r becomes shifted to some field (j) which is bigger than 0o- This just corresponds to the process of eternal self- 
reproduction of the inflationary Universe studied in Section III Dl 



Eq. (|66|l suggests that the total volume of the part of the Universe which experiences one jump to from ^o, grows 
more than exp ^— ^^^2^4 + ci(/io0^ times during the subsequent classical rolling of the field (f) back to 4>q. But this 
factor is much larger than 1 for (p > (p* . Even if such domains later do not jump up but move down according to the 
classical equations of motion, their total volume at the moment when the energy density inside them becomes equal 
to po ^ 10~^^(7 • cm~^ is exp ^— ^^^2^4 + cicpocj)^ times larger than the volume of those (typical) domains which did 

not experience any jumps up at all. And the voltime of those domains which experienced two, three or more jumps 
up will be even larger! 

This means that almost all physical volume of the Universe in a state with given density (for example, on the 
hypersurface po ~ 10~^^g ■ cm~^) should result from the evolution of those relatively rare but additionally inflated 
regions in which the field 0, over the longest possible time, has been fluctuating about its maximum possible values, 
such that V{(f>) ^ 1. 



C. Expansion of the Universe as a measure of time 



In the previous Sections we used the time measured by synchronized clock of comoving observers as a time coordinate 
t. However, in general relativity one may use different ways of measuring time intervals. For example, instead of 
synchronizing clocks of different observers, we may ask them to measure a local growth of the scale factor of the 
Universe near each of them. Namely, the local value of the scale factor in the inflationary Universe grows as follows: 

a(x,t) = a(x,0) exp iJ(0(x,ti)) dti^ . (67) 

Then one can define a new time coordinate r as a logarithm of the growth of the scale factor, 

T = ln^p^ = f H{P{xM)) dt, . (68) 
a{x,0) Jo 

If one neglects the space and time dependence of H, then r = Ht. However, in a more general case the difference 
between t and r may be significant, and, as we will see in the next Section, the time t often is more convenient. 
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Now let us try to derive equations describing diffusion as a function of the new time r. 
First of all, instead of eq. (|49|l we have now 

rf7 ) - ' (69) 

and the ordinary classical equation of motion of a homogeneous scalar field during inflation acquires the following 
form: 

This gives n = and V = = The final form of the diffusion equation is [s^ 



= 3^ 90 1^^^(^) 90 — + "^^^ ^ ) • ^ ^ 

There are two important advantages of this choice of time over the standard one [3^ Is^ . First of all, if we 
consider evolution of a domain with initial value of the field (p — (po, then the physical wavelength Xp of perturbations 
which are frozen at each given moment of time t is proportional to H~^{(j>). This corresponds to the wavelength 

Ac = H~^((j)) exp (^—J^ H{(j3{x,ti)) dtij in the comoving coordinates. The dependence of on (jj is not very 

important. Indeed, the deviation of 4> from (po becomes significant only on an exponentially large length scale. 
However, the second term significantly changes even for small deviations of (f> from (j)Q. Therefore, when we will 
perform our computer simulations of stochastic processes in the inflationary Universe, we will face a complicated 
problem of adding to each other waves of the field (f) with wavelengths which rapidly change from one point to another 
(see next Section). The only way of doing it which we found is to use the recently developed concept of wavelets, waves 
with a compact support, see e.g. [55j . There is no such problem in the time r, since with this time parametrization 

the integral Jq H{(I){x, ti)) dti is the same for all points with given r. 

Another important advantage of using the time t is an extremely simple transition from the distribution Pc(0, t) 
to Pp(0,r): 

Pp(0,T) = Pc(0,T) . (72) 
Naively, one would expect that the normalized distribution Pp(0, r) is trivially related to Pc(0, r), 

Pp{<j>,T) = P,{(f>,T) . (73) 

However, this is true only if the distribution Pc{4', t) is normalized. It does not make much sense to normalize Pc since 
this probability distribution is not stationary, and decreases in time in all realistic models of inflation. Meanwhile, 
the normalized probability distribution Pp may become stationary. Therefore the normalized probability distribution 
Pp will be proportional to Pc, but the coefficient of proportionality will depend on r. 

The first step towards finding the normalized probability distribution Pp is to write a diffusion equation which 
describes the non-normalized distribution Pp H72|) : 

Despite all advantages of using the time r, we should remember that this is not the time which can be measured by 
usual clock of a local observer. Rather it is a peculiar time which an observer measures by his rulers. In what follows 
we will discuss both time t and time r, but we will mainly concentrate on the evolution in time t. 



IV. COMPUTER SIMULATIONS 



All concepts discussed in this paper are rather unusual. We are used to think about the Universe as if it were an 
expanding spherically symmetric ball of fire. Now we are trying to explore the possibility that the Universe looks 
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uniform only locally, but on a much larger scale it looks like a fractal. The properties of fractals are not very simple, 
but there is some beauty and universality in them, which transcends the beauty and universality of simple spherically 
symmetric objects. The simplest way to grasp the new picture of the Universe would be to make a computer simulation 
of its structure. 

Unfortunately, it is very difficult to make any kind of graphical illustration of physical processes in the inflationary 
Universe, since the Universe is curved, and, moreover, it is curved differently in its different parts. What we made is 
just a small step in this direction. We performed a computer simulation of stochastic evolution of the scalar field (jj 
in one- and two-dimensional inflationary Universes, both in time t and in time t ^ loga(i). 



A. Fluctuations in a one-dimensional Universe 



As a first step, we considered a part of a one-dimensional inflationary Universe, of initial size Hq = H~^{(j)Q) (i.e. 
the size of the horizon) . Then we followed generation of perturbations of the scalar field 4> during expansion of this 
part of the Universe. We did it by dividing the process of expansion into small steps of duration At = u H^^{(f>o), 
where u is some small parameter. At each step we solved the equations of motion for classical field (p at all points 
of the grid, and then added to the result a 'quantum fluctuation' Scj). The idea was to represent such quantum 

fluctuations by waves with random phases and directions, with an amplitude \/ //(^t) ' with the (physical) 

wavelength Xp, which would be equal to H^^{(f)) at each point. Note that the amplitude of a fluctuation is equal to 

for u = 1, H = Hq. This gives correct contribution to < {5(p)'^ >^ after summation over the contributions 

with different phases.^ The extra factor "/^t^^^ appears due to the fact that we add each wave not within the time 

H^^{<t)), but within the time uH~^{(j)o). Introduction of the small phenomenological parameter u provides a better 
model of white noise, since it allows for a possibility of jumps higher than the average amplitude within the time 

(even though the probability of large jumps will be exponentially suppressed). We took u ^ 10^^ in most of the 
calculations. 

All calculations were performed in comoving coordinates, which did not change during the expansion of the Universe. 
In such coordinates, expansion of the Universe results in an exponential shrinking of wavelengths, 

-J' ff(,/.(xi,ti))dti 

Ac(x,t) = \p{x,t) e " . (75) 

Perturbations which have a wavelength in physical coordinates, have the wavelength ^ H^^ in comoving coor- 
dinates. Thus, at each step of calculations we were adding the perturbation 

S(f>(x,t) — a • sin / H((f>(xi,t))--^ — ^— dxi + an 
\Jo Mxi,t) 

/ i?((/)(xi, t)) e " dxi + a„ . (76) 

Here a„ are random numbers, and the amplitude of the perturbation is given by 



H{cf>{x,t)) uHi<l>{x,t)) 

We have made our one-dimensional calculations using the grids containing up to 3 x 10^ points. At each step of 
calculations we added sinusoidal waves H76|) corresponding to perturbations of the field cj) generated during the time 



In fact, the amplitude of fluctuations in the one-dimensional Universe would be slightly smaller l56l . but in our computer simulations 
we preserve its normal (three-dimensional) amplitude. 
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uHf^ ^, and then subtracted the decrease of the field (j) due to its classical motion, 



2 J 2 

In our computer simulations we considered the simplest theory of a massive scalar field with Vicj)) = We 
took the mass m very big (m = 0.5), since otherwise, according to eq. H28|l . we would need exponentially large pages 
to show our results. In most of the calculations the initial value of the scalar field was about 1, which corresponds to 
the energy density which is one order of magnitude smaller than the Planck density. 

The results of the computer simulation arc shown in the set of Figures 1. Some explanations are necessary here. 
We present a sequence of 3 panels corresponding to different moments of time t. Each panel consists of three figures. 
On the first panel we show the spatial distribution of the scalar field (f) in comoving coordinates. The value of the 
scalar field at each point is given by the upper boundary of the shaded area. The distribution of the scalar field in 
these coordinates averaged over the whole domain corresponds to the probability distribution Pc- 

The second panel shows the same distribution of the scalar field, but in different coordinates, which take into 
account different rate of expansion in different domains. To achieve this goal, at each step At of our calculations 
we expand the distance between the nearby points by e^('^(^'*))'^*. After that, the computer displays the resulting 
distribution at the interval of the original size, as if there were no change in the total length of the domain. The 
distribution of the field (p in these figures corresponds to the probability distribution Pp, and the squeezing of the 
interval back to its original size automatically normalizes this distribution. This means that the distribution of the 
scalar field shown at the lower figure gives a good idea of the relative fraction of volume of the Universe filled by the 
field (/). 

Well, this is not quite correct. These figures give a good idea of the relative fraction of length of a one- dimensional 
Universe filled by the field (j). Indeed, our computer simulations here were one-dimensional. To get a correct idea of 
how the probability distribution Pp behaves in a three-dimensional Universe, each time At one should expand the 
distance between the points by g^H(tj}(x,t))\t ^ This gives a distorted spatial distribution of the scalar field, but this 
distribution shows in a correct way the distribution of the three-dimensional volume of the Universe filled by the field 
(f). This is shown at the third panel (the third series of figures in Fig. 1, the lowest ones). 

Black and white regions in these figures correspond to the distribution of another scalar field, which we added to 
our model. This is a scalar field $ with the effective potential V{^) which may have several different minima. It is 
assumed that during the stage of inflation driven by the field 0, the effective potential V{<^) is much smaller than 
V{4)). Therefore, at this stage one may neglect the contribution of the field $ to the rate of expansion of the Universe. 
However, this field may be responsible for the symmetry breaking in the theory of elementary particles. The field 
also experiences quantum fluctuations and Brownian motion. 

This Brownian motion can push the field $ from one minimum of 1^($) to another. Then, after the end of infiation, 
the field $ becomes trapped by the minimum to which it jumped, and it cannot move from it anymore. However, in 
different exponentially large domains of the Universe this field may be trapped in different minima. As a result, the 
Universe becomes divided into exponentially large domains with different types of elementary particle physics inside 
each of them. 

To give a particular example of such theory, we will consider a model with the effective potential 

!/($)- Uo (^1-cos— j , (79) 

where N is some integer. The ratio ^ is considered as an angular variable. Such potentials appear, e.g., in the theory 
of the axion field, as well as in the model of 'natural inflation' driven by the fleld The amplitude of the quantum 
jumps of the fleld is described by the same equation (j76|l as for the fleld 4>] to study the classical motion of the fleld 
$ one should apply eq. (jTSfl to this field. 

In our series of computer simulations of the one-dimensional inflationary Universe, we take a particular value N = 2. 
In this case the effective potential (|751l has two minima of equal depth, V{n^o) = U(— 7r$o) = 0. We use black color 
to show the parts of the Universe with $ ^ tt^o, we use white color to show the parts with <I> ^ —tt^q, and we use 
various shades of grey to show intermediate regions. 

Our calculations begin with a domain of initial size H-^ filled by some homogeneous (or almost homogeneous) 
fields (f) and $. After few steps the distribution of the scalar field (j) becomes slightly inhomogeneous. At the same 
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FIG. 1: Time evolution of the distribution of the infiaton scalar field </>, shown by the upper boundary of the shaded area, 

and the scalar field $, shown by the color (black - white). This distribution is shown as a function of one coordinate x, in a 
domain of initial size H^''{(j)o). We present a sequence of three panels for each time t. The first one shows the distribution of 
the fields in the comoving coordinates, which is related to the probability distribution Pc. The second one exhibits the same 
fields in the coordinates which show the physical distance from one point to another, divided by the total distance between the 
two sides of the domain. The third one shows the distribution of the fields per unit physical three-dimensional volume, which 
would correspond to the distribution Pp in a three-dimensional Universe. 

time the inhomogeneities generated in the field $ may be relatively large, if their amplitude 5$ ~ y/uH is large in 
comparison with $o- A-S a result, the field $ rapidly jumps from one minimum of its effective potential to another, 
and the Universe becomes divided into approximately equal number of black and white regions. 

After few steps of calculation the scalar field (/) also becomes very inhomogcncous. However, this inhomogeneity 
looks different on the three different series of figures we produced. In comoving coordinates (upper figures) the field 
^ in the most part of space decreases; we see only few hills, which in the course of time grow and produce many 
thin spikes. However, in the physical coordinates, which take into account exponential expansion of the Universe, 
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these thin spikes look hke large mountains, see figures in the middle. Gradually these mountains occupy more and 
more space, and their physical size grows exponentially. Also, these mountains are never sharp; they are build from 
the sinusoidal waves which always look very smooth on the scale H^^. (The two last statements are only partially 
illustrated by the figures, since the mountains look more narrow than they really are after the computer shrinks the 
expanded interval back to its original size.) 

Finally, the last series of figures shows a very rapid growth of the relative fraction of volume of the Universe occupied 
by the growing field (p. This exactly corresponds to the process of self-reproduction of the Universe. We see also 
that quantum fluctuations of the scalar field $ lead to formation of exponentially large number of black and white 
domains, corresponding to different types of symmetry breaking in the theory (|79|l . 



B. One-dimensional Universe in the r-parametrization of time 



Similar calculations can be done for the r-parametrization of time, where r = In "l^^gj . This immediately reveals 
an important difference between this parametrization of time and the standard one. First of all, now all three sets of 
figures look absolutely the same, since, by definition, at a given time r the degree of exponential expansion is the same 
for all points x. An important difference between the corresponding distributions Pc and Pp is that the distribution 
Pc is not stationary: the comoving volume of the regions filled by the large field 4> decreases. Meanwhile, the total 
volume of the domains filled by the large field 0, which is given by Pp, increases exponentially. As we will show in 
this paper, after we divide Pp by the overall growth of the volume of the Universe, we obtain a stationary normalized 
distribution Pp. 

Unfortunately, this difference between Pc and Pp cannot be seen in our figures. Therefore we will present only 
one figure (Fig. 2) instead of the three series of figures. This figure shows the distribution of the scalar fields (f> and 
$ after a few steps of our computer simulations. It illustrates a specific difference between the distributions at a 
given t and at a given r. The typical wavelength of inhomogeneities of the field (f) in the r-parametrization of time is 
approximately constant in different parts of the Universe, whereas the field cf) in the t-parametrization is much more 
inhomogeneous near the maxima of its distribution in comoving coordinates. 




FIG. 2: The distribution of the fields (j) and $ after several steps in time r. 

From the technical point of view, the difference between the two calculations is refiected in the equations for the 
perturbations Scj) and 5^, and for their classical shift during the time interval At — u corresponding to one step of 
our calculations. For example, for the field </> we have 



50(x, y) = . sm{H{(p{x)) e^x + a„ 

V 27r ^ 



and 



A0(a::,y) - 



uV'{(i)) 
iH^{<P{x,t)) 



(80) 



(81) 
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C. Fluctuations in a two-dimensional Universe 

When generalizing our methods for two- or three-dimensional Universe in time t, one meets several complications. 
First of all, it is possible to plot the distribution of the field in comoving coordinates, but it is very difficult to stretch 
these coordinates to show the distribution of the fields per unit physical volume, as we did for the one-dimensional 
Universe. This seems to be a rather general problem, since it is not always possible to show the curvature of space 
by representing it as a curved surface in a flat space of higher dimension. Therefore we represented all our results in 
the comoving coordinates only. 

Another problem is related to the possibility to represent quantum fluctuations of the scalar field by plane waves. 
In one-dimensional Universe this was an easy problem, but in a two dimensional curved Universe this is impossible. 
In some sense, the waves of the scalar field look as if they were propagating in a medium with an exponentially large 
refraction coefRcient. This refraction coefficient changes from point to point. If one starts with a plane wave at some 
point X, it does not remain a plane wave at some distance from this point: the front of the wave becomes strongly 
bent because of gravitational lensing. 

This is not a real problem if one studies only relatively small domains, but even if one starts with a very small 
domain, eventually its size becomes exponentially large, and the effects of bending become more and more significant. 
Therefore, if one wishes to investigate evolution of the scalar field in the whole domain, one should find some other way 
to simulate quantum fiuctuations of the scalar field, which would lead to the same (or almost the same) correlation 
functions (|f 9|I - H22|) as the ordinary plane wave fluctuations, while being free from the problems mentioned above. 
Such a method does exist: one should use functions 5(j){x) which look like plane waves inside domains of a sufficiently 
large radius rg, where rg ^ H^^. However, these functions should rapidly decrease at r > tq. Such fmictions are 
called wavelets |55l |. 

A simplest example of a wavelet which we used in our calculations looks as follows: 

50(:E,i) =oexp(-C2(r-r„)'-if'(r„,i))) X sin f i/(r„, t))^4^^(2;sin0„ -f y cos0„ -f a„) ) • (82) 

Here a = \J ^^H{lo f^' ^" ^"^^ ^""^ random numbers, C is some constant, (r — r„)^ is the square of the 

physical distance from the center of the wavelet rn to the point r. In the comoving coordinates x and y 

(r - r„)^ = [{X - x^f + [y - y^f) ■ . (83) 

At C ^ 1 this field configuration looks like a sinusoidal wave with an amplitude which is given by a at r = r„, and 
which becomes exponentially small at a distance about {CH)~^ from the center of the wavelet. In the limit C ^ 
wavelets would behave as ordinary plane waves. Thus, one should take C small, but not too small, to avoid the effect 
of gravitational lensing inside the wavelets. 

The procedure of simulation of quantum perturbations of the scalar field consisted of generation of a proper number 
density of wavelets and of their distribution at random points rn in such a way that the number of the wavelets in 

a given domain is proportional to its physical volume (i.e. to g/ ''*). By choosing a proper number density 
of wavelets we mean that the variance of perturbations produced by wavelets in the limit when the Universe is 
homogeneous should coincide with the standard result < (50)^ >= H^t/Aii^. 

Our choice of the form of the wavelets is not perfect, but it is good enough for our purposes. Ideally, one should 
consider wavelets which form a complete orthonormal set of functions, like it is done in the usual Fourier analysis. It 
is well known that the wavelet transform can give a much better information about very inhomogeneous structures 
than the standard Fourier transform . It is also known that the usual momentum representation, which is related 
to the Fourier transform, is not suitable for formulation of quantum field theory in curved space. We were forced 
to use wavelets instead of plane waves to perform a consistent computer simulation of quantum fiuctuations in the 
usual t-parametrization of time. This gives us a hint that the wavelet transform, rather than the standard Fourier 
transform, may become a very useful tool for field quantization in curved space. 

On the other hand, in some simple but important cases one may escape from the difficulties associated with 
gravitational lensing by choosing a different set of coordinates. In our case everything becomes very simple in the 
T-parametrization of time, where r ~ loga(i). Computer simulations of the evolution in time r are much easier jS^. 
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In particular, one can use ordinary plane waves, but with the amplitude proportional to — ■ The reason is 

that different parts of inflationary Universe expand by the same number of e-foldings at a given value of the time 
r. From the point of view of computer simulations this means that expansion of the Universe does not introduce 
any distortions into the original grid of comoving coordinates x and y.^ However, it would be incorrect to make 
computer simulations in time r instead of simulations in time t. These two approaches are complimentary to each 
other, emphasizing different features of the same process. For example, as we have already seen in our study of 
one-dimensional Universe, sometimes it is easier to illustrate the process of self-reproduction of inflationary domains 
using ordinary t-parametrization of time. Figures 3-6 contained in this section were obtained by computer simulations 
of the evolution of scalar fields in time t, whereas Fig. 7 corresponds to the time r. 




FIG. 3: Tlie evolution of the scalar field <jf> in a two-dimensional inflationary Universe in time t. 

The set of Figs. 3 shows the evolution of the scalar field in time t. We begin our simulations with an almost 
homogeneous scalar field (f> ^ (f>o inside a domain of initial size 7J~^(0o)i which is represented by a grid containing 
1000 X 1000 points. Fig. 3.1. Then the Universe expands, the amplitude of the classical scalar field slowly decreases, 
but due to overlapping of waves (=wavelets) corresponding to fluctuations of the scalar field, the amplitude of the 



^ There still remains a small distortion, since the physical wavelength ~ H~^{<)>) takes different values in different parts of the Universe. 
However, this is a relatively small effect as compared with the exponential distortion discussed above; in computer simulations one may 
simply keep the wavelength proportional to H~^{ij>o). 
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scalar field in some parts of the expanding domain becomes much greater than cjjo- This corresponds to growing 
mountains on Figs. 3.2-3.6. We should emphasize again that the peaks of these mountains are not sharp at all; they 
are built from fluctuations with the physical wavelength ^ H^^ . They look sharp since we are plotting everything 
in comoving coordinates, in which the wavelength of perturbations Ac becomes exponentially small, especially in the 
parts of the Universe where <j) is large and expansion is fast (i.e. on the tops of the mountains). 

In the parts of the Universe where the scalar field becomes small there are no high mountains; the field continues 
moving towards (/) = 0. We live in one of such parts of the Universe. The perturbations of the scalar field generated 
at this stage of inflation are relatively small. These perturbations are responsible for the small perturbations of 
temperature of the microwave background radiation discovered by COBE. 

However, if 0o is greater than some critical value (/>*, then mountains with (p > (jf are permanently produced. This 
corresponds to the process of the Universe self-reproduction. The effective potential of the field </) on tops of some of 
these mountains becomes greater than the Planck density. (In these simulations we did not introduce any boundary 
conditions at the Planck density; see however the next Section.) One may envisage each of such tops as a beginning of 
a new "Big Bang" . If one wishes to reserve this name for the first "Big Bang" (if there was one), one may think about 
such names as a "Small Bang" or "Pretty Big Bang" . Whatever the words, the theory we have now is considerably 
different from the old Big Bang theory. 




FIG. 4: Generation of axion domain walls during inflation. 



The series of Figures 4 corresponds to a more complicated theory of elementary particles, which includes another 
scalar field, $, with the effective potential (|79|) with 1. In this case the shape of the effective potential V{^) 

coincides with the shape of the effective potential in the simplest version of the axion theory. In these figures we 
do not show the evolution of the inflaton scalar field 0; thus, no mountains. We just show the distribution of the 
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fluctuating scalar field $ on the plane {x, y). White regions correspond to the minima of V{^) at $/<l>o = 2n7r, black 
regions correspond to the maxima at $/$o — (2n- + I)'''. 

In the beginning of our calculations the color (the level of grey) only slowly varies in the domain of initial size 
~ Fig. 4.1. Later on, in the main part of the domain (in comoving coordinates) the field $ moves towards the 

minima of its effective potential, but in some places the field diffuses, say, from the minimum at $/$o = to the 
minimum at $/$o = 27r. On the way it passes through the maximum at $/$o = ti", and in some places the field stays 
at this maximum, forming a domain wall between the domains with $/$o = and with <i>/$o = 27r. The existence 
of such domain walls in the axion theory was first pointed out by Sikivie |57| . The inflationary mechanism of their 
production described above was discussed in The series of Figures 4 shows the process of formation of these 

domain walls, which look as thin black lines at Fig. 4. 




FIG. 5: Axion domain walls in a model with a smaller value of $o (Pollock Universe). 

The distribution and the properties of domain walls in the axion theory are very sensitive to the values of parameters 
of the theory. For example, by reducing the radius (f>o of the effective potential (|79() one obtains much more domain 
walls per unit volume, see Fig. 5. For obvious reasons, we called this figure 'Pollock Universe'. 

Now let us consider the same model, but with TV = 2, as we did in our one-dimensional calculations. In this case 
in the interval < ^/^o < 27r we have two different minima with V — 0. From the point of view of perturbations 
produced during infiation, there is no much difference between the cases iV = 1 and N = 2; domain structure will 
be generated in each case. But now we assume that some properties of particles interacting with the field $ change 
when the field <&/$o jumps from one minimum to another. This happens, for example, in the supersymmetric SU{5) 
model, where there exist many minima of the same depth corresponding to different types of symmetry breaking in 
the theory. To distinguish between the minima at <I>/<I>o — 2mT and the minima at <I>/<I>o = (2n + l)7r we will paint 
the first ones white and the second ones black, just as we did in our one-dimensional calculations. 

The results of the calculations are represented by a series of Figures 6. Originally the whole initial domain contains 
the scalar field in one of the minima of its effective potentials. Then this field begins jumping from one minimum to 
another, and the Universe becomes divided into many exponentially large black and white domains. Figures 6 show 
that in the course of time the domain which originally was quite homogeneous becomes looking like a huge fractal. 
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FIG. 6: Generation of domain structure of the Universe during inflation in the axion model with N = 2 (Kandinsky Universe). 

We called this fractal 'Kandinsky Universe'. 

In Figs. 7, which we obtained by studying the of stochastic processes in the time r, we plot simultaneously the 
distribution of the inflaton field and of the field By looking at these figures one can find that the most chaotic and 
rapidly changing distribution of black and white domains corresponds to the regions near the peaks of the distribution 
of the field (j). In these regions the color of the domains changes very rapidly, which means that the laws of the low- 
energy elementary particle physics (which are related to the type of spontaneous symmetry breaking) are not fixed 
there yet. However, in the valleys, where the inflaton field (p is small, the color becomes fixed. We live in one of such 
domains of a given color. Other domains are exponentially far away from us. And there are domain walls separating 
domains with different colors. Laws of low energy physics change as one jumps through the domain walls; one should 
think twice before doing so. 

An interesting feature of all these computer simulations is the following. If one takes a magnifying glass and looks 
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FIG. 7: The distribution of the scalar fields (j> and $ in a two-dimensional Universe after several steps in time r. 



at some part of the picture of a physical size of the order of H^^, one will see a rather homogeneous distribution 
of the scalar fields (j) and <&. Then one may start computer simulations again. And again one will get mountains of 
the same type, and a Universe consisting of domains of different colors. This is the first (and the most fundamental) 
kind of stationarity which we certainly have in our Universe. Not only the Universe looks like a fractal at any given 
moment of time, it actually is a growing fractal, which reproduces itself over and over again. 

As we already mentioned, the fractal structure of the inflationary Universe in the theories where inflation occurs 
near a local maximum of effective potential V{(j)) is relatively simple; one can easily calculate the fractal dimension of 
the inflationary Universe js^]. The corresponding structure in the chaotic inflation scenario with V{4)) ~ or with 
V{4)) ~ e°"^ is much more complicated, but it is also much more interesting. From the point of view of computer 
simulations, the fractal structure of the Universe in this scenario has more "colors" in it, since quantum fluctuations 
near the Planck density allow more profound changes of the vacuum state. 
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V. TOWARDS THE THEORY OF A STATIONARY UNIVERSE 



A. Are there Any Stationary Solutions for Pp? 



Now let us return to the problem of stationarity of the distribution Pp{(j),t). Eq. suggests that there are 

no stationary distributions not only for Pc, but for Pp as well. This conclusion, however, may fail when the field 
approaches (pp and density approaches the Planck density pp ~ Mp = 1. One of the reasons is that at p > 1 the 
probability distribution Pc becomes strongly non-Gaussian, dispersion of fluctuations of the field (j) becomes very large 
and our derivation of eq. (|56|l does not work. To avoid this problem, one may try to find exact solutions of eq. Ht)5|l . 

An attempt to do this was made in very important papers by Nambu and Sasaki |3S| . It was claimed that in 
theories V{(j)) ~ with n < A the distribution Pp is always nonstationary, but in theories with steeper effective 
potentials {V{4)) ^ 0" with n > 4, or ¥{({)) ^ e"^"^) stationary distributions for Pp do exist. However, the stationary 
distributions which were found in |38| have their maxima at densities p ^ pp — 1. These results are unreliable for 
several different reasons: 



1. Eq. H55|) may have a slightly different form, which corresponds to the difference between possible definitions 
of stochastic term in the equation (|48|l . This difference is not important at densities much smaller than 1, but 
at p ^ 1 it may become significant. (Note, however, that this is not an unsolvable problem. One may just 
investigate modified equations as well. Two other problems are more fundamental.) 

2. Diffusion equations were derived in the semiclassical approximation which breaks down near the Planck energy 
density. 

3. Interpretation of the processes described by these equations is based on the notion of classical fields in a classical 
space-time, which is not applicable at densities greater than 1 because of large fluctuations of metric at such 
densities. In particular, our interpretation of Pc and Pp as of probabilities to find classical field at a given 
point (or in a given volume) at a given time does not make much sense at p > 1. 



Thus, the results of ref. do not help us to establish the existence of stationary solutions. However, it is obvious 
that there exists a class of theories where Pp is stationary even though V{(t)) grows indefinitely at large </>. For example, 

one can consider a theory with the effective potential ~ A0"exp^^^j , with C 3> 1. (Similar effective potentials 

often appear in supergravity.) In such a theory we have both inflation and self-reproduction of the Universe. However, 
at > C(j)* inflation is impossible since the potential is too steep. Therefore the distribution Pp will be unable to 
move to (j) > C(f)* , and will become stationary. 

There is another, more general reason to expect the existence of stationary solutions. As we will argue now, inflation 
tends to kill itself as the energy density approaches the Planck density. 

In our previous investigation we assumed that the vacuum energy density is given by V^cj)) and the energy- 
momentum tensor is given by V{4>)g^,^. However, quantum fluctuations of the scalar field give the contribution 
to the average value of the energy momentum tensor, which does not depend on mass (for <C H^) and is given by 

m 

<T^u>^ Hi .9m- = \v^ 9^.u ■ (84) 



One of the sources of this contribution is obvious. Quantum fluctuations of the scalar field freeze out with the 
amplitude ^ and the wavelength ~ H~^. Thus, they lead to the gradient energy density (9^(50)^ ^ H^. 

Note, that eq. (|84|l does not give the total contribution of quantum fluctuations to the energy density. When the 
field 4> is outside the minimum of the effective potential, additional terms appear, in particular, the term ~ V'{(f>)d(j) 
(|23|l which is responsible for galaxy formation. At < when the density perturbations responsible for galaxy 
formation were produced, the vacuum energy renormalization (|84|l is subdominant. On the other hand, a,t 4> > (p* 
the contribution corresponding to eq. becomes greater than V'{(j))S(j) ~ V'{(j))H/2'K. For example, in the theory 
V = m202/2 
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for > (/)* ^ 1 / -y/m. 

An interesting property of eq. H84|l is that the average value of the energy-momentum tensor of quantum fluctuations 
does not look like an energy-momentum tensor corresponding to the gradients of a sinusoidal wave. It looks rather 
like a renormalization of the vacuum energy- momentum tensor (it is proportional to g^u). This result appears after 
averaging over all short-wavelength fluctuations and over all possible outcomes of the process of generation of long- 
wavelength perturbations. Something similar happens when one calculates the energy-momentum tensor in ordinary 
flat space-time: the contribution of each wave is not proportional to g^j^, but the energy-momentum tensor becomes 
proportional to g^i, after integration over all such contributions with the invariant measure ~ 6{k^ — m?). However, 
if we are not integrating over all long-wavelength perturbations, but treat them as a classical inhomogeneous scalar 
fleld, then at each particular time interval ~ H^^ in each particular /i- region the energy momentum tensor of these 
perturbations is mainly due to the gradient energy density (5^(50)^, and it is not proportional to g^^. 

This does not lead to any interesting effects at y <C 1 (0 ^ 0p) since in this case <C V . However, at the 
density comparable with the Planck density the situation becomes much more complicated. At 1^ > 1 the gradient 
energy density ^ becomes greater than the potential energy density V{4>). A typical wavelength of perturbations 
giving the main contribution to the gradient energy is given by the size of the horizon, / ^ H^^. This means that 
the inflationary Universe at the Planck density becomes divided into many domains of the size of the horizon, density 
contrast between these domains being of the order of one. These domains evolve as separate mini-Universes with 
the energy density dominated not by the potential energy density but by the energy density of gradients of the field 
4>. Such domains drop out from the process of exponential expansion. Some of them may reenter this process later, 
but many of them collapse into black holes within the typical time and then evaporate. Indeed, the standard 
criterion for the formation of primordial black holes is exactly the requirement that the local increase of density 5p 
is comparable to p on the scale of the horizon. This criterion is satisfied for perturbations of density produced at 
^(0) > 1- 

Of course, one may argue that all our considerations do not make sense at densities greater than the Planck density. 
When the energy density in any /i-region approaches the Planck density, it may no longer be described in terms of 
classical space-time and should be just thrown away from our consideration. In particular, its volume should not 
be considered as contributing to the total volume of the Universe. Thus, such domains should be neglected in our 
definition of Pp. In this case the distribution Pp for the field will stop growing and will approach a stationary regime 
when this distribution will be shifted towards 4> 4>p. 

We are making an even stronger statement. Even if one makes an attempt to consider the domains with F > 1 as 
a part of classical space-time, many such domains drop out from the process of inflation. This means that the total 
volume of inflationary /i-regions cannot grow as fast as e^^^*^-** for </> > (pp. 

In the models describing many different scalar fields inffation may become self-destructing even at an energy density 
somewhat smaller than the Planck density. Let us consider, for example, an inflationary model describing N different 
scalar fields, with ^ 1. The Hubble constant in such a theory is determined by the sum of all effective potentials, 
V — V1 + V2 + .... The Planck density can be defined as a density at which all higher order gravitational corrections are 
equally important. The one- loop correction to the energy- momentum tensor (|84|l in our model is directly proportional 
to the number of scalar fields but the second-loop corrections do not contain such an enhancement. Thus, one 
may argue that the Planck density in our model remains the same, pp ^ 1. On the other hand, the energy density of 
gradients of the scalar field fluctuations now is proportional to N: 

id^6(j))^ NV\ (86) 

This means that the gradient energy density becomes greater than the potential energy density V at V ^ 1. 
In the model with large N this does not automatically imply black hole formation, since density perturbations will 
be suppressed by I/^/N. However, if the energy density of any domain of a size / ^ becomes dominated by 

the gradient energy, such a domain instead of inflation enters the regime of a slow power-law expansion. At this 
stage, previously generated fluctuations of the scalar flelds either oscillate or at least considerably decrease. This 
breaks down the standard scenario of the new scalar field fluctuations freezing out on the top of the previously frozen 
fluctuations. 

In fact, some predecessors of self-destruction of inflation show up already at sub-Planckian densities even at small 
N . When we discussed the derivation of the diffusion equation for Pc, we assumed that all perturbations of the 
wavelength give rise to inflationary domains. However, this is not the case for very large fluctuations of the 
field (j) with the amplitude of a fluctuation 50 > 1. Such fluctuations should be present if, as we assumed, the term 
^(t) in the Langevin equation (|48|l corresponds to the white noise. The gradient energy of the perturbations of the 
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wavelength H ^ with the ampHtude 50 > 1 is given by ^ (df^Scj))^ > ~ V{(p). Domains where such jumps occur 
may drop out of the process of exponential expansion. Suppression of inflation in the domains of a size produced 
by large jumps effectively implies that the effective noise in the Langevin equation is not entirely white. ^ 

The probability of such a large jump within the time H^^ is suppressed by exp (^^^^^j^r—j ■ This indicates that at 

H 1 (i.e. at sub-Planckian densities) the cut-off of noise corresponding to fluctuations with the amplitude 6(1) > 1 
is not very important, since the probability of such fluctuations is exponentially suppressed. However, the closer we 
are to the Planck density, the stronger is the deviation of £^{t) from the white noise, the less efficient is the stochastic 
process producing new inflationary domains. 

What are the possible consequences of this effect? First of all, we know that if the initial value of the field <j) is in 
the interval (j)* < 4> < 4>p, then the distribution Pp((j),t) moves towards greater and greater 0. The only reason why 
it happens, despite the fact that the distribution Pc(0, i) moves towards small </>, is an additional exponential growth 
of the number of inflationary /i-regions due to expansion of the Universe. If, as we argued, the number of inflationary 
/i-regions with (p (j)p grows at a much slower rate (or does not grow at all, which is the case if we just do not consider 
domains with > 0p as belonging to our classical space-time), then the distribution Pp(0,t) stops moving towards 
large as soon as the maximum of this distribution approaches <j)p. Thus, the distribution Pp{4>,t) may approach 
some stationary regime, being concentrated at sub-Planckian densities. Now let us study this possibility at a more 
quantitative level. 



B. Self-Reproduction of Inflationary Domains as a Branching Diffusion Process 



When describing the process of self-reproduction of inflationary domains one should keep in mind that this is not 
an ordinary diffusion process, but process where diffusion of the scalar field in each particular /i-region is accompanied 
with their branching into many independent /i-regions. As was argued in |59j |. agood candidate for the mathematical 
model describing such behavior is the theory of branching diffusion processes j6(] | . 

There are two main sets of questions which may be asked concerning such processes. First of all, one may be 
interested in the probability Pp(0, t|0o) to find a given field at a given time t under the condition that initial value 
of the field was equal to some 0o — 4>{t = 0). In what follows we will denote 0o as x- 

On the other hand, one may wish to know, what is the probability Pp{(j), t\x) that the given final value of the field 
(j) appeared as a process of diffusion and branching of a domain containing some field x- Or, more generally, what are 
the typical properties of branching Brownian trajectories which end up at a hypersurface of a given 0? 

In the end of the Section IlII Bl we have shown that the probability to find a given field (p in a given volume moves 
towards Planckian densities with the growth of time t, if ipQ > </)*. We have shown also that a typical trajectory 
producing domains of a given density for a long time fluctuates at very large 0, close to the Planck density. Thus, 
we obtained qualitative answers for the two questions mentioned above. To obtain a more quantitative description 
of the process of self-reproduction of the Universe, we will first develop more general methods of investigation using 
the simplest characteristics of the branching diffusion processes which involve only those parameters of the model |53 
that can be unambiguously identified at this stage of our studies. 

Many properties of a branching diffusion process can be reproduced with the help of the following function 

uix,t,V) = {Kt,T^)\ ^o = x) ■ (87) 

Here /i(i,I?) is the number of "particles" with the coordinates x in the interval T) at the moment t. In our case 
the analog of a particle with a coordinate x is the /i-region with a given field 4>. The interval of all possible values of 
(j) corresponding to inflation in classical space-time is Vmax = [0e, Ppl- However, one may choose to consider as V in 
eq. (|87|) any smaller segment of the interval [0e, (pp]. 

By {^{t,'D)\ (po — x) we denoted the mathematical expectation value (or, in other words, the mean value) of the 



A more adequate statement would be that the noise remains white, but when it occasionally becomes "too loud" , it destroys Brownian 
particles. 
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number of "particles" in the interval 2? at the moment t under the condition that initially there was one "particle" 
with (f)Q — X- 



The last condition implies that 

u(x,0,I?) = 



1, ifxeP , 

0, otherwise 



There is some subtlety in definition of the 'number of particles' fi{t, V) in the stochastic approach to inflation. 
Depending on the way one formulates the problem, this may be either the number JV{t, V) of /i-regions or their total 
volume V{t,'D) divided by some normalization factor H~^{(j)^)J The difference between these two quantities stems 
from the fact that the total volume of each /i-region is proportional to H^'^{(p), where ip is the local value of the scalar 

field in a given /i-region. This difference typically is not very important: it is proportional to the third power of "^^y j 
whereas most of the distributions we will obtain depend on H{(j3) exponentially. Nevertheless, whenever appropriate, 
we will make a distinction between these two understandings of V). In what follows, we will understand by /x(t, V) 
the volume V(t, V) of the /i-regions with the scalar field S P. (For an equation describing the number of /i-regions 
Af see also ^|). 

Now we will try to find an equation describing the time evolution of u{x,t,'D). The space-time structure of the 
branching diffusion process may be represented by a growing tree of world lines of diffusing particles. The section 
of the tree at some level represents the spatial distribution of the particles at the corresponding moment of time, 
while the increasing number of "branches" corresponds to the increasing number of particles due to the processes of 
"splitting" of one "parent" particle into several "daughter" ones. Here we produce a derivation of the reduced form 
of the equation describing the growth of this tree |59l | in a way which does not involve any new model parameters in 
addition to those already presenting in the chaotic inflation. 

Consider the whole tree of the branching diffusion process from t = to i + At and divide it into two parts — a 
part from t — to t — At and a part from t = At to t-l- At. After the time At, the original domain flUed by the scalar 
field X — 4>o — = 0) becomes filled by the shghtly inhomogeneous field 0(At), and its volume grows exp(3i?(x)At) 
times. To first order in At, the new domain can be represented as a combination of two domains: the domain of 
original volume V containing the field 0(At), plus the new domain of the volume V • 3-ff(x)At with an unchanged 
field 4> = X- Very soon, after the characteristic time ~ H~^, the interior of the second domain will loose any contact 
with the interior of the first one ("no-hair" theorem for de Sitter space), and its subsequent evolution would proceed 
totally independently. This means that when we take averages calculating u(xit + they split into two parts 

(in the first order in At): 



i{x,t + At,V)^{u{cl){At),t,V)\ cf>o^x)+3H{x)At-u{x,t,V) . (89) 



This yields 



9 . . t i-n^^r Wi^t),t,V)-uix,t,V)\ c^o = x) ,QnA 
— u(x, t, V) = 3H{x) ■ u{x, t, V) + Jim^ — . (90) 

According to the theory of stochastic processes , the last term in this equation can be represented in the following 
way: 

Jim^ i<mt\t.V)-<X.t,V)\^.=x) ^ ^^^^^ ^) ^ ^g^^ 

where the operator A (the generating operator of diffusion) can be constructed in a standard way if the corresponding 
Langevin equation (|48|l is known: 

Af( ^ g'^'(x) d / 3/2. ,df{x)\ V'ix) dfix) 

^^^^^ ^ r " mx)^ • ^''^ 



''' We denote volume by V to distinguish it from the effective potential V. 
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Thus, the function u(x,t,2?) satisfies the following equation 

^ u{x, t,V)^A u(x, t, V) + mix) ■ u(x, t, V) . (93) 

As we have discussed, in our approach this function gives the expectation value of the volume V of all /i-regions 
with the scalar field (j> E V. Now we will consider the interval I?(0) — [<j>e, 0] and define the probability distribution 

PpWAx) = -Q^ ■ (94) 

One can easily understand that this is the same (unnormalized) distribution Pp((j),t) which we studied in the 
previous section, but from now on we will also keep track of the dependence of this function upon 0o = X- In 
particular, according to eqs. this function satisfies equation 

^^PM^t\x) = APp{<P,t\x) + 3H{x) ■ Pp{cl>,t\x) , (95) 

or, in an expanded form. 

This is the branching diffusion analog of the backward Kolmogorov equation known for the ordinary diffusion. The 
probability distribution Pp(0, i|x) satisfies also the forward Kolmogorov equation (i.e. Fokker-Planck equation) 

^ Pp(<^, t\x) = it Pp(0, t\x) + 3H{cP) ■ Pp{(t>, t\x) , (97) 

where A'^ is the operator which is adjoint to the diffusion generating operator A H92II |60l |. From the standard definition 
of the adjoint operator, 



F(0)i ficf>) dcp^ j {it F{c^)}f{c^) d(t> , (98) 
one can easily obtain the following expression for i^ 

^'/<«^^(^J^(-"^«H-^H- 

Thus, the forward Kolmogorov equation is 

This equation coincides with the equation l|65|) . which was earlier derived in a different way. 

C. Boundary conditions 

We will give now a brief discussion of boundary conditions imposed on Pp{(j), t\x)- A more detailed discussion will 
be given in the forthcoming paper [4^ . 

As we will see, in many models the distribution Pp{4>, t\x) is sharply peaked at a very large field <j} ^ <j}p ^ (pe, and 
the maximal value of Pp{4>^ t\x) is much greater than Pp[(j}e, t\x)- In such models the form of boundary conditions at 
is almost irrelevant, and we may simply assume that 



= 0. (101) 
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A similar condition can be imposed on Pp{(f>,t\x) at x = (f>e. 

However, this simple trick does not have a universal validity for all potentials V{(j)) and all possible time parametriza- 
tions. Therefore we should find a more general way to impose boundary conditions at (f>e- 

Note that at the boundary (j) = inflation ends, there is no diffusion back from the region (j) < (j)e, and the 
field (j) in the domain with (j) < <f>f. continues rolling down to even smaller values of </>. The best way to describe 
it is to say that the diffusion coefficient vanishes at (/> < </)e. This means that the form of the diffusion equation 
changes discontinuously at (j) — ipe- However, neither the probability distribution nor the probability current can be 
discontinuous at this point. 

The simplest way to describe this situation is to consider first the distribution Pc, for which we have the probability 
conservation equations (|53|l and l|54|) . (For a description of this method see |6lj.) 

The continuity condition for the probability distribution Pc and the probability current Jc at 0e can be written as 
follows: 

Pc(</>e+) =Pc(</'eJ; Jc{4>e^) ^ Jci^e^) . (102) 

The notations (f),,^ and 4>e- are used to show that we may approach (jy^ either from (f) > (j)^ or from (j) < <f>e. The last 
condition can be written as 



iV'{<t>) p., X 3^(0) 
81/1/2(0)^=^'^^+^ 8^1/2(0) 



.)■ (103) 



Here we have used our assumption that diffusion does not contribute to Jc at < 0e- For generality we restored the 
parameter (5, corresponding to the ambiguity of the definition of stochastic force. 

Using the continuity condition Pc{4'e+) — Pc{4'e-)j we obtain the following boundary condition for P^. 

d 



or, equivalently. 



, (104) 



The stationary solutions we are going to obtain will not be very sensitive to the choice of /3, unless we choose it 
extremely large. In the particular case /3 = i, which we consider throughout the paper, the boundary condition is 

^c(0ej 3 V 



A similar approach can be used for the case of the probability distribution Pp. Here we have a slight complication, 
since the analog of the probability current in this case is not conserved due to the production of new /i-regions. This 
is precisely the reason for the appearance of the additional term 3HPp in the diffusion equation. However, one can 
derive the same boundary conditions in this case as well. In order to do it, one should just integrate equation 1)64(1 
over (j) in the infinitesimal interval from 0e_ to (j)e^ ■ Since the term SHPp is not singular in this interval, we recover 
eq. pU3|l for Pp, and finally obtain the same boundary condition as before, 

^Pp(0ej=-^y Pp(</'eJ. (107) 

Finally, we should fix the boundary conditions for the backward Kolmogorov equation ((96|l at x = Xe = ^e- 
There exists a regular way to reconstruct these boundary conditions from the boundary conditions for the forward 
Kolmogorov equation (|107|l [6l] |. One should accurately reconstruct the diffusion operator A from its adjoint H99|l . 
taking into account that the functions on which these operators act may not disappear at the boundaries. This will 
produce additional boundary terms in the expression for A. Then one should take such boundary conditions for the 
backward Kolmogorov equation, which would make these additional terms vanish. As we will show in a subsequent 
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publication, the resulting boundary condition has a very interesting and general form, which does not depend either 
on the parameter /? or on a particular choice of time parametrization: 



This equation implies that the probability distribution near Xe is always given by the square of the tunneling wave 
function, 



In fact, as we will see soon, in some models to be considered this remarkable relation holds (in a sense to be discussed 
later) in a very large region, almost up to the Planck boundary x 4'pi ^-nd this result remains valid even if we 
considerably modify the boundary condition (|108|) . 

Now we should discuss the boundary conditions near the Planck boundary, where V{4>p) ^ 1. As we have shown in 
the Section IV Al the situation at this boundary is much more complicated and ambiguous. Description of evolution 
of the scalar field at (/) > 0p in terms of classical space-time is impossible, and our diffusion equations do not make 
much sense there. There are several possibilities to be considered. 

1. First of all, we may exclude from our investigation all /i-regions which jump into the state with <j) > <f>p, since we 
cannot study their evolution in any case. In other words, we may restrict our attention to the part of space which 
at all times t > Q remains classical. There are several formal ways to do so. The simplest one is to introduce 
a phenomenological description of destruction of all domains which jump to the region (j) > <j).p. This trick will 
effectively remove these domains from our consideration. This can be done by the same way we introduced 
branching. To do it, we added the term 3H{(j))Pp to the diffusion equation. Now we may, for example, multiply 
this term by a function F{(j)), which is equal to 1 for (j) < <f)p^ and which becomes large and negative for ip > (pp. 
By doing so, we do not prohibit the diffusion process in the region cf) > (pp, but we discard all domains which 
diffuse there. 

2. There is also another possibility which we discussed in the previous section: inflation may lead to a self- 
destruction of inflationary domains with (p > cpp. Phenomenological description of this regime may be achieved 
by the same trick as in the previous case, except for the behavior of the function F{(p) at p > pp may be less 
dramatic; it may just introduce a small coefhcient in front of 3iJ(0). However, if this coefficient is small enough, 
the final effect will be basically the same: it makes the distribution Pp stationary and concentrated at (p <, (pp. 
This would also correspond to something like an absorbing boundary conditions (or, more generally, to "elastic 
screen" type boundary conditions, which are intermediate between the absorbing and reflecting ones). 

3. Finally, one should take into account that in the realistic theories the effective potential may be steep at large (p. 



We have already discussed one of such possibilities, which appears if the effective potential is ~ X<p" exp(^^=^ j , 

with C > 1. In such a theory we have both inflation and self-reproduction of the Universe. However, at p > C(p* 
inflation is impossible since the potential is too steep. Therefore the distribution Pp will be unable to move to 
p > C(p* and will become stationary. 

There may be other reasons why inflation cannot penetrate into the region p > 1. It well may be so that space-time 
with p > 1 simply cannot exist [6^ |6^ , or that the stringy nature of interactions does not allow us to penetrate 
to distances smaller than the Planck scale [63. l65l|. Whatever the reasons are, to describe their phenomenological 
consequences for the distribution Pp one should either add some terms, corresponding to destruction of inflationary 
domains to the diffusion equation, or to modify the effective potential V , the diffusion coefScient T) and the mobility 
coefficient n at p > (pp, or to impose some kind of absorbing or reflecting boundary conditions so as to prevent 
penetration of the fluctuating flcld into the domain with (p (pp. 

Our investigation has shown 43] that the flnal results for the distribution Pp are not very sensitive to the method 
one uses to prevent penetration of the field (p into the domain with p^ pp and on the type of the boundary conditions 
imposed (whether they are absorbing, reflecting, etc.). They depend only on the value of the field p> = p>b, where the 
boundary conditions are to be imposed, and this dependence is rather trivial. Our arguments suggest that p)b ^ (pp. 




(108) 




(109) 
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In what follows, when solving the diffusion equations we will simply assume that the functions u{x,t,T>{4>)) and 
Ppi4'TAx) satisfy absorbing boundary conditions at (/> = ipp, where V{(j)p) = 1. In particular, we will assume that 

Ppicl^p,t\x) ^ , (110) 

where V{(f>p) = 1. A similar condition should be imposed at the Planck boundary for the field x' 

Pp{cj),t\xp)=0 , (111) 

where Xp = 4'p- However, we will keep in mind for future discussion that, strictly speaking, (j)p in these equations is 
just a phenomenological parameter, not necessarily corresponding to VYc/jp) = 1. A detailed investigation of different 
boundary conditions near the Planck boundary will be contained in |43| . 



D. Stationary Solutions 



There are several ways to proceed further. One may try to obtain solutions of equations l|96(l and l|l()()(l directly. 
We prefer a different strategy, looking for a solution of eqs. (I96f) and H1UU|) in a form of the following series of 
biorthonormal system of cigcnfunctions of the pair of adjoint linear operators A + ?>H and + ?>H: 

oc 
s=l 

Indeed, this gives us a solution of eqs. (|96|l and H10Q() if 



2 2^ dx\ 277 dx ' 7 3i?(x)rfx 

and 

The orthonormality condition reads 



r<t>p 



i^s{x)'^Ax)dx = Ssj . (115) 



In our case (with regular boundary conditions) one can easily show that the spectrum of Xj is discrete and bounded 
from above. Therefore the asymptotic solution for Pp{(j), t\x) (in the limit t — > oo) is given by 

Pp(0,t|x) =e'^*^i(x)7ri(0) • (l + 0(e-(^i-^^)*)) . (116) 

Here tpi{x) is the only positive eigenfunction of eq. (|113|l . Ai is the corresponding (real) eigenvalue, and 7ri((/)) 
(invariant density of branching diffusion) is the eigenfunction of the conjugate operator pi4|) with the same eigenvalue 
Ai. Note, that Ai is the largest eigenvalue. Re (Ai — A2) > 0. This is the reason why the asymptotic equation H116() 
is valid at large t. We found that in realistic theories of inflation a typical time of relaxing to the asymptotic 
regime. At ~ (Ai — A2)~^, is extremely small. It is only about a few thousands Planck times, i.e. about 10"^" sec. 

Now we see that the average volume of the Universe filled by the inflaton field (f> in the interval V = [(pe, 4'p], grows 
exponentially. And the exponent e^^* does not depend on and x- This means, that the normalized distribution 

Pp((/.,i|x) = e-^i*Pp((/.,i|x) (117) 
rapidly converges to the time-independent distribution 



Pp{(j>,t^ oo|x) = -01 (x) ((A) 



(118) 



40 



It is this stationary distribution that we were looking for. The remaining problem is to find the functions ■01 (x) ^-nd 
TTi ((}()), and to check that all assumptions about the boundary conditions which we made on the way to eq. are 
actually satisfied. The boundary conditions on Pp in terms of tti and can be written as follows: 

^liXe) = I ^ ■ MXe) , MXp)^0, (119) 

and 

(0e) = ~ J y • ^1 ('/'e) , TTi (</)p) = . (120) 



The general program of finding solutions for V'l and tti for a wide class of inflationary models will be pursued in 
^. In what follows we will present solutions of these equations for the theories V = jcj)^ and V = Vq e°"^. 

Equations (|113|1 . (|114|l for V'i(x) and ni{(j)) in the theory jtj)'^ look as follows: 



,// ,// 6 3\ , /367r Ai /67r\3/2\ , „ , 



, / 6 9\ / 6 15 367r Ai /67r\3/2N 



AA05 0/ VA06 02 A04 TT^nA 
The corresponding boundary conditions pi9|) and p20|) for this theory are 



^i(Xe) = 1^0l(Xe) , V^l(Xp)-0, (123) 

and 

^l(0e) = -|-7ri(0e) , ^l(</)p) = 0, (124) 

/ \ 1/4 

where we take for definiteness (pe ^ Xe = 0.3 and 'Pp — Xp — { j) ■ (Note the difference between the coupling 
constant A and the eigenvalue Ai.) 



The analytic solution of this equation shows that in the limit of small A the eigenvalue Ai = 2v67r w 8.681 |43| . 
if one identifies the upper boundary (f>p with the value of (p at which V{4i) — 1. However, this limit is approached 
very slowly. We solved this equation numerically with the boundary conditions H123fl . H124|l . We have found the 
eigenvalues Ai corresponding to different values of the coupling constant A: 



A 


1 


10-1 


10-2 


10-3 


10--* 


10-^ 


10-*^ 


Ai 


2.813 


4.418 


5.543 


6.405 


7.057 


7.538 


7.885 



As a test of self-consistency of our approach, one may solve one of the equations first, find the eigenvalue Ai, and 
then check that the same eigenvalue can be obtained from the second equation with the boundary conditions specified 
above. We performed this check, and confirmed our results. 

One can find also the second eigenvalue A2. For example, for A = IQ-^ one gets A2 — 6.789. This means that for 
A — 10"* the time of relaxation to the stationary distribution is At ~ (Ai — X2)^^ ^ 4Af-i ~ 10"*^ seconds — a 
very short time indeed. One should note, however, that the complete time for establishing the stationary distribution 
depends on initial conditions, and in some cases it may be much longer than (Ai — X2)^^ 43] . 

Note that the parameter Ai shows the rate of exponential expansion of the volume filled by a given field (f). This 
rate does not depend on the field cp, and has the same order of magnitude as the rate of expansion at the Planck 
density. Indeed, Ai should be compared to 3H{(j)) = 2\/67rF(0), which is equal to 2-\/67r at the Planck density. As we 
already mentioned, in the limit A the eigenvalue Ai also becomes equal to 2^/6Tr « 8.681. The meaning of this 
result is very simple: in the limit A ^ our solution becomes completely concentrated near the Planck boundary, 
and Ai becomes equal to ZH{<pp). 
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At first glance, independence of the rate of expansion of volume e'^^* on the value of the field (f) may seem coun- 
terintuitive. The meaning of this result is that the domain filled with the field (j) gives the largest contribution to 
the growing volume of the Universe if it first diffuses towards the Planckian densities, spends there as long time as 
possible expanding with nearly Planckian rate, and then diffuses back to its original value (f). 

But what about the field (f> which is already at the Planck boundary? Why do the corresponding domains not grow 
exactly with the Planckian Hubble constant H{(j)p) = 2\/67r/3? It happens partially due to diffusion and slow rolling 
of the field towards smaller (p. However, the leading effect is the destructive diffusion towards the space-time foam with 
(j) > (j)p. One may visualize this process by painting white all domains with V{(l)) < 1, and by painting black domains 
filled by space-time foam with V^(0) > 1. Then each time H^^{(/)p) the volume of white domains with ~ 0p grows 
approximately e'' times, but some 'black holes' appear in these domains, and, as a result, the total volume of white 
domains increases only e^'^'i/^^/^ times. This suggests (by analogy with calling the factor djr — 3Ai/2-\/67r 'the 
fractal dimension of classical space-time', or 'the fractal dimension of the inflationary Universe'. (Note that dfr < i 
for A 7^ 0; for example, dfr = 2.6 for A — 10^^.) However, one should keep in mind that the fractal structure of the 
inflationary Universe in the chaotic inflation scenario in general is more complicated than in the new or old inflation 
and cannot be completely specified just by one fractal dimension j43j . 

The numerical solution for 7ri((/)) is shown in Fig. 8. It has some interesting properties. As we already mentioned, it 
is concentrated heavily at the highest allowed values of the inflaton field. This concentration becomes more and more 
pronounced with a decrease of the coupling constant A. At small values of the field (p this function rapidly vanishes. 
This result is in a contrast with the behavior of the square of the Hartle-Hawking wave function, which is extremely 
sharply peaked at small (j). Note that the shape of the solution, as well as the eigenvalues Ai, prove to be practically 
independent on the boundary conditions at = 0e- 



^i(x) 




FIG. 8: Functions TTi{(f>) and 'tpiix)^ representing the stationary solution Pp{cj>,t\x) for the theory j(j>'^ for A — 0.1. 

The function tpiix) looks similar to tti (</>), see Fig. 8, but its functional dependence on x is, in fact, quite different, 

revealing an important similarity to the square of the tunneling wave function ~ exp(^—3/8V{x)j- This is not 

unexpected: the boundary conditions H1U8|I . 1)1091) suggest precisely that kind of behavior. However, this 'easy 
explanation' is, to some extent, misleading. The solution for '01 (x) is not very sensitive to the boundary conditions 
X — Xe, and its relation to the tunneling wave function has deeper reasons, to be explained in [i^ . 

To verify (and clarify) this statement by numerical methods, we will consider the function ^'i(x) related to the 
function Vi as follows: V'i(x) = *i(x) ex.p(^-3/&V{x)j = 'I'i(x) cxp^^-^^)- According to lfnT|) . ifT^ . this 
function obeys the equations 



VAx'^ 



/367r 



12 
A^ 



Ai /67r\3/2 



, 



(125) 



*'i(Xe)=0, 'fi{Xp)=Q. 



(126) 



The solution is shown in Fig. 9. For definiteness, we normalized ^'i(xe) = 1- One can easily see that, e.g., for 
A = 0.3, the function ^'i grows ~ 10^^ times before it reaches its maximum. However, in the same interval, the 
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FIG. 9: Behavior of the function 



*i(x)=^/'i(x) exp(3/8F(x) 



for the theory 



for A = 0.3; time t. 



function exp^—3/8T^(x)^ = exp^— ,^"^^4 ^ grows ~ 

10i°'^2 tijnes. This means that the leading contribution to V'l (x) in 

the whole interval from Xe to Xp is given not by ^Pi, but by the square of the tunneling wave function. This amazing 
fact is in a complete agreement with a tentative interpretation of the tunneling wave function: Its square gives the 
probability that a given domain of an inflationary Universe originally (at t = 0) was (created) in a state with a given 
field (t){t = 0) = x- 

However, one should not be too excited about it. The most dramatic growth of tpi occurs at relatively small 
values of x, close to Xe- In this region the square of the tunneling wave function correctly describes V'l- Similarly, 
the change of ipi in the interval from Xe to Xp in the leading approximation is given by the square of the tunneling 
wave function. This means that the square of the tunneling wave function reasonably well describes the exponential 
suppression of probability that the inflationary Universe was originated in a state with a very small field X ^ Xe- 
However, ed, x ~ Xe ^ Xe the tunneling wave function changes slower than \E'i(x). In particular, the tunneling wave 
function becomes almost constant near the maximum of "01 ■ Therefore if one wishes to find the behavior of V'l near 
its maximum (which is one of the most interesting problems), the tunneling wave function becomes irrelevant, and 
instead of it one should find a complete solution for (x) ■ This was one of the main purposes of our investigation. 



Now let us find the distribution Pp{<j),t\x) for the theory with the exponential potential V{(j>) = Vq 







corresponding equations in this case are 

'3a 3a 



8Vn 



97r 



~\2Vo 



3/2 



-3qx/2 



The 
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9a 3a 



-a0 



9a' 



3a' 



9tt 



-arp 



4 8V0 J V 8 16Vb Vb 
The boundary conditions (|119(l . H120|) in this case look as follows: 



Ai / Stt 



3/2 



-3a0/2 



= 



(128) 



V'KXe) = — V'l(Xe) , V'l(Xp) = , 



(129) 



and 



3a 



7ri(0e) , 7ri((/>p) = 



(130) 



Without any loss of generality, one may take Xe 



in these equations, see next paragraph. The solutions of 



these equations look qualitatively similar to the solutions to equations for the theory see Fig. 10. However, here 
one should issue a warning, which is not very important with this time parametrization, but will be more important 
in the next Section. 

The boundary conditions at the end of inflation in the theory ^j- were natural in the sense that inflation by itself 
ended at (/> = 0e. In the theory Vb e"* inflation never ends by itself; one must change the shape of the effective 
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^i(x) 




FIG. 10: Functions Tri{(j>) and V'i(x)i representing the stationary solution Pp{(j>,t\x) for the theory Vb e"**" with Vb — 0.2, 
a = 0.5. 



potential 'by hand' at some place, say, at = 0. Without such a cut-off one may get all kinds of unphysical stationary 



solutions. The cut-off may be effectively performed, e.g., by considering models with the potential 4Vo sinh^ This 
potential at large positive looks like Vb e°"^, but at small (p it looks like ^rn^cj)^ (where = 2a^Vo), which implies 
that inflation ends at ip <^ 0.2. 

Alternatively, one may introduce an abrupt cut-off a,t (j) — (pe- This can be accomplished by a discontinuous increase 
of the derivative V'{<j)) at (p < (f>e- However, such a change would lead to an extra term in the boundary conditions 



(ITIlHll : 



P. 



8yi/2(^) 



(l^'(</'eJ-^'(0e+))Pp(0e+) 



(131) 



Since the solutions we found extremely rapidly decrease at small 4>^ this modification was not important for us. 
However, with different time parametrizations, the corresponding solutions may take greater values at 0e or Xe- In 
such case one should either take into account improved boundary conditions, or consider potentials of the type of 
4Vb sinh^ ^ instead of the exponential potential Vq e"*^. We will return to this question in the next Section. 



VI. STATIONARY SOLUTIONS WITH OTHER TIME PARAMETRIZATIONS 



In this section we will briefly describe our results concerning the stationary regime in the r-parametrization of time, 
where r = In ^^^^ ~ /J 7?(0(x, ti)) dti, see Section llH CI In this case equations analogous to (|113|l and H114|) look 
as follows (compare with eq. I|74|l 'l: 



and 



1 d 
Stt d(b 



TTjiP)) 



(3-A,)7r,(0) 



(132) 



(133) 



As before, the asymptotic solution for Pp{(j), t\x) (in the limit t — > oo) is given by 

Pp{c^, r\x) = e^'^ Mx) M4>) ■ (l + O (e-(^i-^^)-)) 



(134) 



Here ipiix) is the only positive eigenfunction of eq. (|132|l . Ai is the corresponding eigenvalue, and 7ri((/)) is the 
eigenfunction of the conjugate operator H133|) with the same eigenvalue Ai. The stationary normalized solution we 
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are looking for is given by 



Pp{<f>,T^ oo|x) = Vi(x)'ri((?i) 



(135) 



where Pp((/>, t|x) = Pp(0, t\x) e-^^\ 

Interestingly enough, the boundary conditions for V'i(x) in the new time parametrization remain unchanged. How- 
ever, due to the change of the expression for the probabiUty current in the r-parametrization, the boundary condition 
on 7ri((^e) now looks slightly different: 



= 0, 



(136) 



or, equivalently, 



1 V 



(137) 



Let us show how the equations for Vi(x) and 7ri(^) look in the theory V{(f)) = j(f>'^. To simplify these equations, 
we will make a change of variables, ^ = A^/'*x, ip = and use the notation ai = 12'k\~^/^ (3 — Ai). Then the 

corresponding equations acquire the following form: 



6 



(138) 



The boundary conditions are 



and 



V'i(ee) = 7F^lfe) . V'l(ep) = 



7^1 (Ve) = 7I"i((?ie) , 7ri((/>p)=0. 



(139) 



(140) 



(141) 



The Planck boundary in the new variables corresponds to = tpp = V2. Inflation in the theory ^(^^ ends at (pg 0.3, 
which corresponds to = = 0.3A^/^. 




0.03 0.0300001 0.0300001 0.0300001 



FIG. 11: Behavior of the function 7ri((p) near ipe in the theory j<j>'^; time r. Here ifie = X^^'^(t>e, A = 10 
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We begin with the equation for tti. As in the previous Section, its solutions prove to be rather stable with respect 
to the boundary conditions at ip^. This is illustrated by Fig. 11. In the beginning the curve goes down, in accordance 
with the boundary conditions H14HI . However, almost immediately it turns up and approaches the asymptotic regime 



(142) 



where C, c are some normalization constants. The change of the slope of this curve at by a factor of two would 
just slightly modify the value of ip where the curve reaches its asymptote (I142II . 




FIG. 12: Functions 7ri(<^) and 'i/'i(C)) representing the stationary solution Pp[ 

if = X^l-^(t) and i = A^/^x- 



, t|x) for the theory ^(^^ in terms of variables 



The result of the numerical solution of eq. 1)1391) is shown in Fig. 12. Interestingly enough, the asymptotic solution 
(|142|l gives an excellent approximation (with an accuracy of few percent) to the exact solution for T^i{p) in the whole 
interval from (fe to ipp. We have found that the parameter ai corresponding to the eigenvalue Ai is given by 41.95. 
This leads to the following expression both for the eigenvalue Ai and for the 'fractal dimension' dfr (which in this 
case refers both to the Planck boundary at (pp and to the end of inflation at (pe)'. 



Ai — dfr — 3 



41.95 
127r 



\/A« 3- 1.1 VA . 



(143) 



The solution of eq. H138|l for "01(0 is also shown in Fig. 12. Note that due to our redefinition of variables X ^ 
— > all solutions here remain the same for all values of A. To study this solution in a more detailed way, it is very 

instructive to represent it in the form ■01 (0 = ^i(C) 6xp|^— 3/8y(^)^ = ^i(C) exp^— 2|r^, as we did in the previous 

Section. According to H138|) . H140|) . this function obeys equations 
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(144) 
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The solution is shown in Fig. 13. As before, we normalized 5'i(fe) — 1- One can easily see that, e.g., for 
A = 0.3, the function 5*1 (C) grows only 7 times before it reaches its maximum. In the same interval, the function 

exp^— 3/8V^(^)^ = cxp^— grows ^ 10^°^^ times. For A — 10^'*, the function 5*1 grows only 2300 times, whereas 

the function cxp^— 3/8V^(^)^ grows ^ times. This difference is much more dramatic than in the standard 

t-parametrization of time. It reflects the fact that in almost the whole interval from to f ~ 0.65, where the function 
01 (^) reaches its maximum, it is correctly described by the square of the tunneling wave function. In terms of the 
original variable x? this implies that in about one-third of the whole interval from Xe = 0.3 to Xp = (4/A)"'^/^ the 
function ipi (x) is given by the square of the tunneling wave function: 

3 \ / 3 



^i(x)-exp(-^)=exp(-^) 



2Ax^ 



(146) 
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FIG. 13: Behavior of the function = i/>i(^) exp^3/8V(C)^ for the theory f 0'' for A = 0.3; time r. 



Thus, the square of the tunnehng wave function does play an extremely important role in quantum cosmology. 
A more accurate fit to the function which can be used in the whole interval from Xe to Xpi is given by 

(~8T%j) (v(x)+o.4 ^ o)- 



exp 



The resulting stationary distribution is 
Pp{(j),T ^ oo|x) ~ exp 
= exp 



(-84))(FWTa4-n)^*"''K-'P-^'>*1 
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--i^) .0exp(-3.5VA</>^ 



(147) 



This expression is valid in the whole interval from (j)e to 0p and it correctly describes asymptotic behavior of Pp{4>, t\x) 
both at X ~ Xe and at x ~ Xp- 

A similar investigation can be carried out for the theory V{(t)) — Vq e°"^. The corresponding equations are 



, fa 3a _^^\ ^ 3^(3- Ai) ^_,^ 

Vn 



(148) 
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The (naive) boundary conditions are 



3Q,g-aXe 

tP'liXe) = — JTT ^l(X£ 
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) , i^i(xp) = , 



and 



n[{(l)e) = - - TTi{(f>e) , ni{(f>p) = . 



The solution for -01 (x) in the whole interval from Xe to Xp looks as follows 
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The solution for 7ri((/)) in the whole range from to (f)p with an accuracy of few percent is given by the following 
simple function: 



(153) 
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The solutions are shown in Fig. 14. Note that now the solution for 7ri(0) is concentrated near the end of inflation, 
at (/) ~ 0e. The reason why we obtained such a solution is related to the 'prescribed' nature of the end of inflation 
in the theory with the exponential potential: Motion of the field is slow and diffusion is large until the very end of 
inflation a,t <j> — <j)e = 0- If the speed of motion of the field cj) near the end of inflation had been vanishingly small, the 
distribution tti (0) would be equal to the square of the Hartle-Hawking wave function. For the exponential potential, 
however, this is not the case, see eq. H153(l . 

As we mentioned at the end of the previous Section, when one describes inflation in the theory with the exponential 
potential, some care should be taken of boundary conditions. One should either modify the boundary conditions 
for the case of the exponential potential, or modify the effective potential itself. For example, one may study the 
theory with the effective potential 4Vb sinh^ However, this potential is considerably different from the exponential 
potential in a very large region < < a^^. This leads to a profound modification of the solutions, shifting the 
maximum of 7ri((/i) towards ip > S> 1. On the other hand, there is another way to modify the theory, by adding 
to it one more scalar field, which triggers a phase transition with an instantaneous end of inflation at </) = 0. Such 
a theory is described in [23 • In this case the potential of the field cj) remains exponential at > 0. In order to 
obtain a phenomenological description of the end of inflation in this model one may assume that the derivative V' 
abruptly increases when the field (jj becomes negative. In this case the boundary condition for the function tti becomes 
(compare with (|131|l 'l 



+ 



where in our case (p^ ~ 0. The r.h.s. of this equation gives a positive contribution to the value of ^[{(j)^) in eq. (|151|l . 
The value of this contribution depends on the jump of V' at the point (f>e. Numerical solution of the corresponding 
equation for different values of this jump shows that the solution remains unchanged in the main part of the interval 
from 0e to (pp, except for a small vicinity of the point (f),,, see Fig. 14. This again confirms that our solutions are 
rather robust with respect to the change of the boundary conditions at the end of inflation. 

A complete stationary probability distribution is given by 



= K"'" O'pi-^p) (e-«-V«) ■ {e-'*-V,)e-'*''' . (155) 



This expression gives a rather good approximation for Pp{4>^ r — > oo|x) for all (p and Xj except for a small vicinity of 
the point (p,, mentioned above. 



Note, that the distributions which we obtained do depend on the choice of the time parametrization. This is not 
unexpected; each event may look differently being described in different coordinate systems. However, it would be 
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very desirable to find an invariant description of our results. Fortunately, the most important qualitative results 
we discussed, the existence of the self-reproduction of the Universe and the existence of a stationary probability 
distribution Pp, do not depend on our choice of t- or r-parametrization of time. 



VII. COMMENTS AND INTERPRETATION 

The results obtained in this paper may seem strange and sometimes even counterintuitive. Therefore we are going 
to discuss their interpretation and an alternative derivation in a subsequent publication [4^ . However, it is necessary 
to make some comments right now. 



A. Do we need Planck density to have stationarity? 

There are two possible kinds of the stationarity. First of all, self-reproduction of inflationary domains implies that 
even in a very distant future there will be many inflationary domains in the Universe containing all possible values of 
scalar fields compatible with inflation. According to the no-hair theorem for de Sitter space, each such domain of a 
radius greater than H^^ (/i-region) will expand practically independently of the processes in the nearby domains. This 
means that the Universe will repeatedly reproduce inflationary domains, which will have statistically same properties 
as the similar domains produced billions of years ago. 

This kind of stationarity is the most fundamental. Its existence is related only to the existence of the regime 
of self-reproduction. In all models we considered in this paper self-reproduction of inflationary domains occurs at 
(f) > (j)*, where V{(1>*) ^ 1 |^. Self-reproduction occurs also in the models where inflation is possible when the field (j) 
is near a local maximum of its effective potential, at V{(j)) 1 |37l l40l | : such models originally were used in the new 
inflationary Universe scenario. Thus, the very existence of the process of self-reproduction of inflationary domains, 
and, consequently, the existence of the first, most fundamental stationarity, does not depend on unknown processes 
at V{(f>) > 1. We will call this stationarity local, or micro stationarity, to distinguish it from the global stationarity, or 
macro stationarity, which refers to stationarity of probability distributions over the whole Universe. 

This second kind of stationarity may or may not exist. Its existence does not follow from any general considerations. 
This is similar to the situation in quantum statistics, where the description in terms of a macrocanonical ensemble 
sometimes is impossible despite the existence of a good description in terms of a microcanonical ensemble. On 
the other hand, in the situations where the Universe is globally stationary, the description of its evolution can be 
considerably simplified. 

As we have seen, in all realistic inflationary models the probability distribution in comoving coordinates Pc is not 
stationary. Fortunately, the probability distribution Pp is stationary in many interesting cases. 

First of all, it is stationary in the models where the inflaton field is fluctuating near a local maximum of its effective 
potential 37, 40]. This is already extremely important. However, for the reason to be discussed in this Section, 
it would be most interesting to obtain stationary solutions in the theories where inflation is possible very close to 

V{4>) ~ 1. 

In this paper (see also have found stationary solutions for Pp in a class of models of chaotic inflation with 

the effective potentials |-(/)^ and Vq e"*; a generalization for other theories will be considered in a separate publication 
|43j |. The main assumption which we made to obtain these solutions is that the self-reproduction of inflationary 
domains is impossible (or at least is strongly hampered) at the density higher than the Planck density. We gave three 
different reasons why this assumption may be reasonable: 

1. Diffusion equations for Pp and our interpretation of their solutions it terms of the distribution of a classical field 
in a classical space do not work for V{<l)) > 1. Therefore it seems that at the present level of our understanding 
of the Planckian physics the only thing one can do is to study only infiationary domains with T^(0) ^ 1 and 
discard those domains which jump to V{(j)) ^ 1. 

2. Large energy density concentrated in the spatial gradients of the scalar field fluctuations hampers the process 
of self-reproduction of inflationary domains with U(0) > 1. 
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3. The only natural mass scale near the Planck density is the Planck mass. Therefore even if the effective potential 
is not very steep at small energy density, nothing can protect it from becoming very curved due to quantum 
gravity effects near the Planck density. In such a case inflation (or the process of self-reproduction of inflationary 
domains) will cease to exist at V{(j)) ^ 1. 

These arguments suggest that inflation (or at least self-reproduction of inflationary domains) cannot exist at V{(j)) > 
1. Therefore one should impose some boundary conditions which do not allow Pp to penetrate deeply into the regions 
with V{(j)) > 1. This immediately leads to the existence of stationary solutions for Pp concentrated at V{(j)) ^ 1. 
The exact form of these solutions does depend on the processes near the Planck boundary. However, as we will show 
in a separate publication |43j |. in many cases this dependence is rather trivial and does not change the qualitative 
behavior of the stationary distribution Pp. For example, one can show that for each (j)p there exists such (/)p ~ (j)p 
that imposing absorbing boundary conditions at gives the same solution for tti (0) as imposing reflecting boundary 
conditions at 0p. This means that the choice between absorbing and reflecting boundary conditions is equivalent to 
the corresponding redefinition of the position of the Planck boundary. Thus, even though a complete understanding 
of physical processes near the Planck boundary is important for obtaining exact stationary solutions for Pp, we do 
not expect that the qualitative features of these solutions in the region (f> <ti (j)p will be dramatically different from 
those which we obtained in the present paper. 

One should keep in mind that in some theories the effective potential may become very steep a,t (p > <f>b where 
V{(j)h) ^ 1, i-e. without any relation to quantum gravity effects and the Planck boundary. If this happens at c/fh > (jj* , 
then we will have a stationary distribution which is concentrated at (j) ^ (jfb, and all our results will remain qualitatively 
correct after substituting (f>b instead of (j)p. In particular, there will be a fractal geometry with a dimension djr = Xi 
(for T-parametrization of time), where Ai will be slightly less than 3 (and, correspondingly, dfr = h{4)^) ^ ^ ^"-"^ 
t-parametrization) . 

On the other hand, there are some theories where the effective potential at large (j} approaches some constant value 
Vo <C 1. In such models nothing will prevent the distribution Pp from moving towards indefinitely large without ever 
reaching the Planck boundary, and the global stationarity may be absent. Another class of models which may have 
runaway solutions are the models including interaction terms ^(jpR. In such models the distribution Pp is concentrated 
near the Planck boundary, but the Planck boundary is a line rather than a point, and there is a room for a runaway 
diffusion along this boundary |66|. 



B. Is it possible to have stationarity and self-reproduction of the Universe at < 0*? 

In the first part of the paper we have shown that self-reproduction of the Universe occurs only if we have an 
inflationary domain with the field (j) > (jf . The critical field </>* is typically much greater than 0e- For example, 
(/)* l/yjrn 3> 0e in the theory ^(jp'. However, we obtained a stationary probability distribution for all in the 
interval from (j),, to (pp. Does this mean that now we are taking our words back and saying that the self-reproduction 
of inflationary domains may occur at <C 0* as well? 

The answer to this question consists of two parts. First of all, for the existence of stationary distribution Pp{4>, t\x) 
at some field (j) — po one does not need inflation and self-reproduction of the Universe at pQ. The distribution will 
remain stationary even for <C 0ej i-e. after the end of inflation. The stationarity of distribution at small (j) is an 
automatic consequence of the stationarity of distribution at large (j) and of the diffusion and rolling of the inflaton 
field from large towards small 0. This is the reason why we obtained the same speed of growth of volume e^^* for 
all 0. This implies also that the total volume of all parts of the Universe with p ~ 10~^^g-cm^^ should increase with 
the same speed. Indeed, according to our results, the total volume of the domains with = 0e at all times increases 
as e"*"^*. But 10^° years later, the density of matter inside all such domains will become ~ lO^^^g-cm""^. It is obvious 
that since the total volume of the domains with = 0e permanently grows as e^^*, the total volume of the all domains 
with p ~ 10~^^g-cm~'^ will grow with exactly the same speed (with the time delay of 10^" years). Consequently, the 
relative fraction of volume of the parts of the Universe with any given properties (i.e. the total volume of the parts 
with the given properties, divided by e^^*) is time-independent, even if these parts are post-inflationary. 

Independently of this question, it is interesting to address the issue of self-reproduction of the parts of the Universe 
with < 0*. 

The critical field 0* plays an important role when we study the possibility that the process of self-reproduction of 
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inflationary domains begins from one domain of a size 0{H~^). Indeed, if we have one /i-region with the field (p > (j)* , 
this is already enough to ensure the permanent self-reproduction of inflationary domains. However, soon after the 
beginning of this process there will be exponentially many domains with all values of the field cj), including the field 
(f> ~ (f>e- Typical quantum jumps of the field (j) in such domains within the time At ^ have the amplitude ^ 
which is much smaller than the average classical decrease A(j) of the scalar field during this time 

V V 

A(j> = = . (156) 

The end of inflation in a typical scenario is dcflncd by the condition that the decrease of V{(p) within the Hubble time 
becomes comparable with V{(t)): |At/((/>)| V'\A(j)\ - V{(j)). This gives V ^ V&t^V, and |A0| l/^/87^ at the end 
of inflation. In all realistic models of inflation this decrease is much greater than the typical amplitude of fluctuations 
at the end of inflation, tP- <C 1. 

However, if we already have many domains with <f> ^ (pf,, then in some of these domains large jumps with an 
amplitude Stj) ^ \A<j)\ ^ I/a/Stt may occur. The probability of such jumps is exponentially suppressed, but once they 
occur and bring the fleld towards ></)*, an inflnite process of self-reproduction of inflationary domains begins again. 

Let us first check whether such jumps are possible at all. There may exist quantum fluctuations with any amplitude. 
However, if the amplitude is too large, then the gradient energy density of these fluctuations becomes much greater 
than the potential energy density V{(j)), and the standard approach we used in this paper should be considerably 
modified j42j. The boundary at which our 'white noise' becomes 'not so white' is determined by the condition 
^{dM? ~ UHS^)^ ^ V{(P), which gives S(t> ^ V6/VS^. 



This means that the standard description of the fluctuations with Stp ^ Acj) is valid for cf) > (pe, since \A<j)\ < I/vStt 
at (/) > 0e. Therefore the portion of the original volume where the fleld cj) within the typical time At — experiences 
a jump up by S<j) = C|A0| ^ C/VStt is given by the Gaussian distribution 

Here C is some constant, which should be somewhat greater than 1 if we wish that the field cj) increases despite its 
decrease due to rolling down by \A(j)\ ~ I/^/Stt: 5(j) > |A0|. If we want the field 4> to experience a subsequent 

jump by the same value Scf) — C/\/8n, we should multiply our result by exp ^" 32V(0+A0+(50) ) • However, each of 
the subsequent jumps will be much more probable than the first one, since the degree of suppression is exponentially 
sensitive to the increase of V{(p + A(j) + Scj)). Moreover, now the field should not jump so high to compete with the 
classical rolling, since the value of A(j) decreases at large (p. Therefore, to a reasonable approximation, the probability 
of the first jump H157() gives us the whole result. (This estimate is particularly good for small initial (p ^ <pe, which 
is comparable to l/^/Sn.) Thus, the portion of the original volume occupied by some field (f){t = 0) = Xi which 
experiences a series of jumps to (p* (or to any other field ip > x), is given by 

/ 3C^ \ 

P{X -^<P)^ exp -— • (158) 



32V{x) 

With an accuracy of the factor C = 2, this is just the square of the tunneling wave function! 

This result implies that if we have an inflationary domain (or a collection of inflationary domains) of a total volume 
~ H~^{(p) exp (^+ 32^v(y) ^ ; some parts of this domain will enter eternal process of the self-reproduction of the Universe 
even if the field inside this domain initially was smaller than (p* . 



C. What about the Big Bang? 



The main result of our work is that under certain conditions the properties of our Universe can be described by a 
time-independent probability distribution, which we have found for theories with polynomial and exponential effective 
potentials. A lot of work still has to be done to verify this conclusion. However, once this result is taken seriously, 
one should consider its interpretation and rather unusual implications. 
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When making cosmological observations, we study our part of the Universe and find that in this part inflation 
ended about te ~ 10"'^° years ago. The standard assumption of the first models of inflation was that the total duration 
of the inflationary stage was At ^ 10~^^ seconds. Thus one could come to an obvious conclusion that our part of the 
Universe was created in the Big Bang, at the time + At ^ 10^° years ago. However, in our scenario the answer is 
somewhat different. 

Let us consider an inflationary domain which gave rise to the self-reproduction of new inflationary domains. As 
we argued in Section 5.2, one can visualize self-reproduction of inflationary domains as a branching diffusion process. 
During this process, the first inflationary domain of initial radius ^ H~^{(j>) within the time H~^{<f>) splits into ~ 20 
independent inflationary domains of similar size. Each of them contains a slightly different field (f), modified both by 
classical motion down to the minimum of V{(j)) and by long-wavelength quantum fluctuations of amplitude ~ H/2t:. 
After the next time step H~^{(l)), which will be slightly different for each of these domains, they split again, and so 
on. The whole process now looks like a branching tree growing from the first (root) domain. The radius of each 
branch is given by H~^; the total volume of all domains at any given time t corresponds to the 'cross-section' of all 
branches of the tree at that time, and is proportional to the number of branches. This volume rapidly grows, but 
when calculating it, one should take into account that those branches, in which the field becomes greater than (j)p, die 
and fall down from the tree, and each branch in which the field becomes smaller than 0e, ends on an 'apple' (a part 
of the Universe where inflation ended and life became possible). 

One of our results is that even after we discard at each given moment the dead branches and the branches with 
apples at their ends, the total volume of live (inflationary) domains will continue growing exponentially, as e^^*. What 
is even more interesting, we have found that very soon the portion of branches with given properties (with given values 
of scalar flelds, etc.) becomes time-independent. Thus, by observing any flnite part of a tree at any given time t one 
cannot tell how old the tree is. 

To give the most dramatic representation of our conclusions, let us see where most of the apples grow. This can 
be done simply by integrating e^^' from t = to t = T and taking the limit as T ^ oo. The result obviously 
diverges at large T as Aj~^ C^^"^, which means that most apples grow at an indefinitely large distance from the root. 
In other words, if we ask what is the total duration of inflation which produced a typical apple, the answer is that it 
is indeflnitely long. 

This conclusion may seem very strange. Indeed, if one takes a typical point in the root domain, one can show that 
inflation at this point ends within a flnite time At ~ 10~^^ seconds. This is a correct (though model-dependent) 
result which can be conflrmed by stochastic methods, using the distribution Pc{(j), At\x) 9]. How could it happen 
that the duration of inflation was any longer than 10~^^ seconds? 

The answer is related to the choice between Pc and Pp, or between roots and fruits. Typical points in the root 
domain drop out from the process of inflation within 10~^^ seconds. The number of those points which drop out from 
inflation at a much later stage is exponentially suppressed, but they produce the main part of the total volume of the 
Universe. Note that the length of each particular branch continued back in time may well be flnite However, there 
is no overall upper limit to the length of branches, and, as we have seen, the longest branches produce almost all parts 
of the Universe with properties similar to the properties of the part where we live now. Since by local observations 
we can tell nothing about our distance in time from the root domain, our probabilistic arguments suggest that the 
root domain is, perhaps, indeflnitely far away from us. Moreover, nothing in our part of the Universe depends on the 
distance from the root domain, and, consequently, on the distance from the Big Bang. 

Thus, inflation solves many problems of the Big Bang theory and ensures that this theory provides an excellent 
description of the local structure of the Universe. However, after making all kinds of improvements of this theory, we 
are now winding up with a model of a stationary Universe, in which the notion of the Big Bang loses its dominant 
position, being removed to the indefinite past. 

But from infiation it also follows that on a much larger scale the Universe is extremely mhomogeneous. In some 
parts of the Universe the energy density p is now of the order of one (in Planck units), which is 123 orders of magnitude 
higher than the ~ 10~^^g ■ cm~^ we can see nearby. In such a scenario there is no reason to assume that the Universe 
was initially homogeneous and that all of its causally disconnected parts started their expansions simultaneously. 

There is one subtle point in this discussion. According to our picture, the main part of the volume of the Universe 
is being produced by inflation not very far from the Planck boundary. Indeed, imagine that we imposed our absorbing 
or reflecting boundary conditions not at V{(t)p) = 1, but at V{(j)p) = 1/4, thus excluding from our consideration those 
branches where the energy density may become greater than V{(j)p) = 1/4. This would reduce the eigenvalue Ai 
approximately by a factor of 2: Ai ~ ^Ai. Correspondingly, the total volume of the Universe would grow not as e'*'^*. 
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but at a much slower rate e'^'^*. This means that the main contribution to the total volume of the Universe is given by 
the branches which from time to time come very close to the state of the maximal possible energy density compatible 
with the self-reproduction of inflationary domains. In Section 4 we called such states 'Small Bangs'; they correspond 
to the peaks of the mountains on Figs. 3.6 and 7.5. 

Therefore even though the original Big Bang singularity may be removed to the indefinite past, the typical time 
lapse from the last Small Bang to the end of inflation may be as small as 10"'^^ seconds. From the point of view of a 
present observer, there is no much difference between the last Small Bang and the original Big Bang; thus one may 
still use the old Big Bang theory for a phenomenological description of the observational data. However, if one wishes 
to understand the global structure of the Universe, its beginning and its fate, one should use a more complete theory 
including quantum cosmology. 

The difference between the Big Bang and the Small Bangs is especially clear in the theories where there exists a 
large difference between the maximal possible value of V{(l)) compatible with inflation and the Planck density. This is 
the case in the theories with the potentials which have a local maximum at some 0, or which approach a plateau with 
V{(j)) = Vo <C 1 at large 0, or which become very steep at </> > 0f,, where <j)* <^ (f>b (f>p. In such theories inflationary 
quantum fluctuations never bring the Universe back to the Planck density. However, as we have argued, this difference 
may exist even in such theories as (/>" and e°"^. The Big Bang is believed to be the single event of creation of the very 
first domains of classical space-time; the description of these domains in terms of classical space-time was impossible 
before the Big Bang. Meanwhile, each branch may come through the state which we call 'Small Bang' many times; the 
more times each branch approaches the state of maximal density compatible with inflation, the greater contribution 
to the total volume of the Universe it gives. 



D. Initial conditions for inflation, different versions of inflationary theory and observational data 

Even though the properties of the Universe around us do not depend on the time from the Big Bang to the present 
epoch, one may still try to examine the problem of initial conditions near the cosmological singularity. However, 
as we argued above, a more relevant problem would be to find typical trajectories (branches) which could produce 
a part of the Universe of our type. Some aspects of this problem can be studied with the help of the probability 
distribution Pp(0, t|x), or with the help of the function V'i(x)- 

Let us consider for example the classical potential of the Coleman- Weinberg type used in the new inflationary 
theory 

VW^X^^Hn^-^ + ^^ . (159) 

Here ^ > 0; the field a <C 1 corresponds to the minimum of the effective potential H159|) . In the new infiationary 
Universe scenario it was assumed that the high temperature effects put the field (p onto the top of the effective 
potential at (j) — 0, and then inflation begins as the temperature drops down Q. However, it was soon realized that 
this scenario usually does not work, and one should consider chaotic inflation instead. According to this scenario, one 
can have two different possibilities to obtain inflation in the theory (|159|l : One may consider inflationary domains 
where the Universe from the very beginning was in a state close to = 0, or domains with > 1 Whereas both 

possibilities in principle can be realized, the second one appears to be much more probable from the point of view 
of initial conditions. As we argued in Section 2, initial conditions which are necessary for inflation naturally appear 
if inflation may begin at V{(p) ^ 1, i.e. close to the Planck density. On the other hand, the probability of infiation 
seems to be exponentially suppressed if it may begin only at V{(f>) <^ 1, which typically is the case if it begins at 
= 0. This conclusion follows from many different arguments, including the qualitative discussion of possible initial 
conditions near the singularity 5, '2Cf|, calculation of the tunneling wave function of the Universe jl3j|-|l6||. and even 
from the computer simulation of different regimes of the Universe expansion |46j |. Did we learn anything new about 
it after the present work? 

First of all, our results give some additional indication of importance of the tunneling wave function, since its square 
appears in the expression for ■01 (x)- This may be considered as a new confirmation of our conjecture that the most 
natural realization of chaotic infiation scenario in the theory (|159|l occurs if infiation begins near the Planck density, 
at A0^1n| - 1, i.e. at > 1. 
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However, now we can say even more. If we examine the solutions of the diffusion equation for Pp in this model we 
would find that these solutions consist of two independent branches. The field cj) in domains with (j) < a never 
diffuses to > fj, and vice versa. Correspondingly, the distribution Pp((/),t|x), which gives the total volume occupied 
by the field 0, will consist of two branches: 

Pp+(^,t|x) =e^^*7r+(^)^+(x) for <^,x>^t, (160) 

and 

p-(0,t|x) =e^i"*7rr(0)Vr(x) for 4>,X<^- (161) 

The eigenvalue is approximately equal to 3H{(f> = 0) ^ 1 H^jEHl, whereas from our results it follows that Af ~ 1. 
This means that if the Universe originally contained both domains with cj) > a and with (f> < a, then within a very 
short time of the order of the Planck time t ^ I, the main part of the volume of the Universe becomes totally 
dominated by domains with (f> > a. 

Using the analogy between the Universe and the growing tree, let us imagine that the branches with (f> > a are 
green, and the branches with (jj < a are red. Our results imply that if the Universe originally contained at least one 
green branch, then very soon it becomes all green, independently of the original number of red branches (which in fact 
should be smaller than the original number of green branches if our previous arguments concerning initial conditions 
are correct). 

This conclusion may be somewhat modified if one takes into account a possibility of large jumps of the field (j) over 
the non-inflationary region near a. Note that the width of this region is very large, A(j) ~ 1, and since this region is 
not inflationary, our methods do not allow any description of such jumps. If such jumps are possible, for example, due 
to the process described in 49], or due to Euclidean tunneling, then green branches will produce red off-springs, and 
the total volume of the red domains will grow with the same speed as the total volume of the green domains. However, 
even if this process is possible, the total volume of domains 4> < a will remain exponentially small as compared with 
the volume of domains with (f> > a because of the exponentially small probability of such large jumps 49|| . 

Suppose now that the effective potential V{(j)) H159|l at large (p grows not as Xcf)* In ^ but, for example, as e°"^ with 
a ^ Vl67r, or in any other way which precludes inflation (or self-reproduction of inflationary domains) at (j) > a. 
Then, in the absence of any competition, the red branches (beginning at <C cr) will win. Of course, this will 
happen only if initial conditions are good enough for inflation to begin at (j) ~ 0. As we just argued, this probability 
is exponentially small. However, we are speaking about the conditional probability that our part of the Universe 
was formed by inflationary expansion, as compared with the probability that it was formed by some non-inflationary 
process. This suggests that in order to compare these probabilities one should multiply the probability of having given 
initial conditions by the subsequent increase of the volume of any part of the Universe with these initial conditions, 
and then normalize the total probability by dividing it by the total volume of the Universe. This was our method of 
obtaining the normalized probability distribution Pp. If this is correct, then the self-reproducing inflationary branches 
always win. Indeed, whatever is the probability of proper initial conditions for inflation beginning at V{4>) <^ 1, the 
main contribution to the total volume of the Universe will be given by self-reproducing inflationary domains, since 
this contribution even now continues growing exponentially as e^^*, and only this contribution will survive after the 
division by e^^*. 

Thus, it may be possible to have a consistent cosmological theory even if inflation occurs only at V{(j)) <^ 1. This 
is the case, e.g., for the 'natural inflation' and for the hyperextended inflation ,51j|. However, as we have seen, 
realization of inflation in these models requires somewhat trickier reasoning than in the theories where inflation is 
possible at V{(t)) ^ 1. What if the Universe was created in a non-inflationary state, because of the exponential 
suppression ^ e-^/^'^W of the probability to produce an inflationary Universe with small V{(l))7 Does it make sense 
to multiply the volume of an unborn inflationary Universe by e"^^*? Even though we believe that this is a wrong 
objection and that the indefinitely large growth of volume produced by inflation does solve the problem of initial 
conditions for every theory where the self-reproduction of the Universe is possible 0, 0| , it would be desirable to 
have an alternative realization of chaotic inflation scenario, where initial conditions for inflation would appear in a 
natural way even if inflation driven by the field cj) may occur only at V{(j)) <^ 1. 

Fortunately, this task can be easily fulfilled [23, . The simplest way to do it is to introduce an additional scalar 
field which may have no interaction with the field (p, but which can drive inflation at V{^) ~ 1. In such a theory 
inflation may begin at V{^) 1 in a quite natural way. This may be a 'bad' inflation, producing large perturbations 
— ^ 10~^. However, in the process of indefinitely long inflation driven by the field <I>, fluctuations of the scalar 
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field (j) are produced, which in some domains of the Universe put the field (p on the top of its effective potential. In 
those domains where the heavy field $ rolls down to the minimum of V{^), inflation continues due to the light field 
0, which determines the properties of the observable part of our Universe. Models of this type were called 'hybrid 
inflation' in ref. [27l| . 

Let us consider, for example, the 'natural inflation' model, which has an inflaton potential of the type 



This is exactly the model which we considered in our computer simulation of production of the axionic domain walls 
during inflation driven by the field (f>, see Section 3.4. Black lines in Fig. 5 correspond to the parts of the Universe 
where $ — 7r(2n+ 1), i.e. where quantum fluctuations have driven the field <i> onto the top of its effective potential. In 
the comoving coordinates these black lines look very thin, but one should remember that they are always much thicker 
than . In those parts of the Universe where the field (f) eventually rolls down to the minimum of its effective 
potential, inflation still continues in the black regions with $ = n{2n + 1). The Universe in these regions enters the 
process of eternal self-reproduction. After a while, the parts of the Universe produced by these black regions enter 
a stationary regime, in which nothing depends on the first stages of inflation driven by the field (j). This exactly 
corresponds to the scenario we discussed in the previous paragraph, up to an obvious redefinition ^ $. 

Of course, there will be many other parts of the Universe where at the end of inflationary stage driven by the field 
0, the fleld $ will not stay on the top of the effective potential (|162|l . The volume of such parts at the end of the first 
stage of inflation will be greater than the volume of the parts with |<i> — 7r(2n + 1)| <^ H hy a. factor ~ <i>o/i/ ^ 1. 
This is reflected by the relatively large area painted white in Fig. 5. However, the regions with |<i> — 7r(2n + 1)| < 
(painted black in Fig. 5) enter the process of self-reproduction, and very soon the main part of the volume of the 
Universe (at a given value of the fleld (/)) will be dominated by domains produced by inflation of these black regions. 

Thus, even though the 'natural inflation' model does not look natural at all from the point of view of initial 
conditions, it can be incorporated into a chaotic inflation scenario of a more general type where the problem of initial 
conditions can be easily solved. According to this scenario, we live in the remnants of the domain walls produced at 
the flrst stage of inflation. 

Our main conclusion is that whenever one may have two different inflationary branches, the main contribution to 
the total volume of the Universe will be given by the branch on which inflation may occur at a greater value of U(0). 
The initial conditions for inflation on this branch are also more natural. However, it is possible to construct consistent 
inflationary models in which our part of the Universe is produced by a branch where inflation may occur only at very 
smah V{(j}). 

The main motivation for our investigation was a desire to get a complete and internally consistent picture of the 
global structure of the Universe. However, sometimes the knowledge of the global structure of the Universe may tell 
us something nontrivial about its local structure. For example, let us consider the ratio of adiabatic perturbations 
to gravitational waves, or the ratio between scalar and tensor perturbations of metric. Is it possible to find a model- 
independent relation between these perturbations |H3,f73|? 

In some models this relation may depend not only on the choice of the potential, but also on the initial conditions 
for inflation. In particular, the amplitude and spectrum of adiabatic perturbations produced at the last stages of 
inflation in the theory (|159|) do not depend strongly on whether inflation begins at <C cr or at </> ^ cr. In both cases 
^ ^ VA log^^^ I 7]. However, the amplitude of gravitational waves produced during the rolling from (f> > a will be 
proportional to i?(0 Mp) ~ 0{V A), whereas the amplitude of gravitational waves produced during the rolling from 
(j) < a is proportional to H{(f> = 0) ~ V^(t>o ^ "v/A- Without saying anything about initial conditions and about the 
process of self-reproduction of the Universe in this theory, one cannot have a definite prediction of the amplitude of 
gravitational waves and of the ratio of tensor perturbations of metric to the scalar ones in this model. However, as 
we have argued, the main part of the volume of the Universe is produced by inflation along the branches with <f> > a. 
This enables us to make a deflnite prediction of the amplitude of tensor perturbations in this model: this amplitude 
should be proportional to \/A rather than to \/A0o- 

This example shows that the investigation of the global structure of inflationary Universe, which at the flrst glance 
deals with the scales which we will never observe, may lead to testable predictions. In what follows we will discuss a 
much more speculative possibility, which, however, may deserve future investigation. 




(162) 
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Until now, most of the studies of properties of the observable part of the Universe were to some extent related to 
the distribution P^. It was assumed that one should take a typical sample of inflationary Universe, study its evolution 
and predict its properties. However, our investigation shows that the field distributions which seem typical from the 
point of view of Pc, may happen to be very unusual from the point of view of Pp, and vice versa. What if we live in 
a part of the Universe which is very unusual from the 'normal' point of view, reflected in the distribution Pc? 

We have discussed already the most dramatic difference between Pc and Pp. A trajectory which is typical from the 
point of view of corresponds to the duration of inflation ~ 10^"^^ seconds, whereas from the point of view of Pp a 
typical trajectory spends an indefinitely long time in the inflationary regime near the Planck boundary. Analogously, 
the standard investigation of density perturbations suggests that the typical deviation of the density of the Universe 
from the critical density cannot exceed ^ ~ 10^^ on the scale of the horizon. The question we would like to address 
is whether there may exist some reason for us to live in a not very typical domain, where the deviation of density 
from the critical density might be greater than 10"^. The reason why this might happen can be explained as follows. 
If one tries to flnd out those trajectories which give the dominant contribution to the total volume of the Universe 
within a given time interval t, one may conclude that if the time t is sufficiently large, then the main contribution will 
be given by the trajectories which first rush up towards the highest possible values of V{4>), spend there as long time 
as possible, and at the very last moment rush down towards small with the speed even somewhat greater than the 
speed of the classical rolling. Indeed, even though the probability of such a regime is exponentially suppressed, those 
lucky trajectories, which happen to follow these rules of the game and saved some time in the process of rolling down, 
may become generously rewarded by the additional exponentially large growth of volume during the extra time they 
spend near the highest possible values of V {(/)). 

If this is the case, the main contribution to the volume of the Universe at a given time t in a state with a given 
density p will be given by the rare regions where the field (j) at the last stages of inflation jumped down with the 
speed exceeding its speed of the nearby domains. Thus, the energy density inside these regions should be somewhat 
smaller than the energy density in their neighborhood. However, not every jump down may compensate the delay in 
the homogeneous classical rolling. The proper jump should result in a nearly homogeneous decrease of density on a 
scale comparable to the scale of the horizon. A region where such a jump took place will look now like a part of an 
open Universe. Such regions will be spherically symmetric, since the probability of formation of asymmetric regions 
by quantum jumps will be suppressed by a large exponential factor, just as the probability of formation of asymmetric 
bubbles during the first order phase transitions*. 

Needless to say, this scenario is extremely speculative. Moreover, a preliminary investigation of this effect indicates 
that in the simple theories with ^(0) ~ </>" the corresponding effect is insignificant. The only theories where this 
effect might be considerable are those where the Hubble constant at the last stages of inflation is very small. However, 
the possibility that we may live in a locally open part of inflationary Universe, i.e. in a part with 51 < 1, should 
not be overlooked. In order to study it in a proper way, one needs to know not only the probability distribution Pp, 
but also the correlation functions for the values of the field (j) at different points, taking into account different rate of 
expansion of the Universe along different inflationary trajectories. We hope to return to a discussion of this possibility 
in a separate publication. 



VIII. BRIEF SUMMARY 



The history of the development of cosmology seems to follow a very nontrivial path. At the beginning of the century 
many people tried to flnd a stationary solution of the Einstein equations, with the hope that General Relativity would 
resolve the inability of Newton's theory to provide us with a stationary cosmological model. Einstein even introduced 
the cosmological constant into his theory for this purpose. The non-stationary character of the Big Bang theory 
advocated by Gamov on the basis of Friedmann cosmological models seemed very unpleasant to many scientists in the 
flfties. Then, the discovery of the cosmic microwave background turned the situation upside down. Physicists began 
to treat with contempt any attempts to flnd stationarity (remember the 'steady-state' model). After several decades 
of the reign of the Big Bang theory, the inflationary scenario appeared, which solved many of the intrinsic problems 
of the Big Bang cosmology and apparently removed the last doubts concerning its validity. 



For a discussion of a related possibility to explain homogeneity of the observable part of the Universe see also 1491 l7ll . 
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However, it was realized soon afterwards, that inflation is even more dynamic than the old Big Bang theory. 
In inflationary cosmology, in addition to the ordinary classical evolution of the Universe governed by the Einstein 
equation, quantum mechanical evolution proves to be extremely important, being responsible for the large-scale 
structure formation and even for the global structure of the Universe. This quantum mechanical evolution can be 
approximately described by stochastic methods, and some of the solutions of the corresponding stochastic equations 
prove to be stationary! Surprisingly enough, after the dramatic development of the Big Bang theory during the last 
ten years, we are coming now to a new formulation of the stationary cosmology, on a new level of understanding 
and without losing a single achievement of our predecessors. The observable part of the Universe can be very well 
described by the homogeneous isotropic Big Bang model. However, on extremely large scales (far beyond the visible 
horizon) the Universe is very inhomogeneous. On even larger scales this inhomogeneity produces a kind of fractal 
structure, repeating itself on larger and larger time and length scales. The statistical properties of this structure are 
what we have found to be stationary. 

Of course, we are far from giving final answers to fundamental questions of cosmology. In fact, we arc just beginning 
to learn how to ask proper questions in the context of the new cosmological paradigm. We need to make sure that our 
results have correct quantum mechanical interpretation. It would be important to find an invariant way of formulating 
our results, which would not depend on a particular choice of time parametrization. There are many other problems 
related to our approach which are to be solved. We know, however, that at the first stage of any investigation in 
quantum field theory or in quantum statistics it makes a lot of sense to find a vacuum state, or a state of thermal 
equilibrium, or any stationary state which may play the role of a ground state of the system. We hope very much 
that the existence of the stationary regime of the evolution of the Universe described in this paper may help us to 
find a proper framework for the future investigation of quantum cosmology. 

It is a pleasure to thank J. Bond, S. Hawking, M. Mijic, D. Page, A. Starobinsky and A. Vilenkin for valuable 
discussions. This work was supported in part by NSF grant PHY-8612280. 
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